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FOREWORD 


This  technical  report  has  been  prepared  by  Dr.  E.  R.  0.  Eokert,  Professor  of 
Mechanical  Engineering  at  the  University  of  Minnesota  under  Contract  AF  33(616) - 
2214  for  the  Aeronautical  Research  Laboratory,  Office  of  Aerospace  Research,  Wright- 
Patterson  Air  Force  Base,  Ohio.  The  work  reported  herein  was  accomplished  on  Task 
70138,  "Investigation  of  Heat  Transfer"  of  Project  7064,  "Research  on  Aerodynamic 
Flow  Fields".  Mr.  Erich  Soehngen,  ARL,  was  the  project  scientist. 

This  technical  report  supersedes  WADC  TR  54-70  dated  April  1954. 

Appreciation  is  extended  to  Miss  Carolyn  Byers  for  typing  the  manuscript. 


ABSHACT 

this  report  la  a  roiaaoo  of  NIDC  It  5i*-7D  ahloh  praaantad  a  aarvoy  of  tha 
literotara  on  baat  tranafar  at  high  onperacnic  velocities  pabllahad  predooimtely 
la  the  yoara  W5l  to  19$Jt.  Tha  larga  nandiar  of  papara  aada  it  naoeaeary  to 
rostriot  tha  aumy  ooeantlally  to  mrfaoes  along  aid  oh  tha  praaaura  la  oonatant, 
aa  aaMpttoi  ahloh  la  not  too  rostrlotiva  alaea  evidence  ooold  ba  presented  that 
for  daartar  ahapes  aa  gantrally  aaad  la  Md  apaad  aircraft  tha  Inflaanaa  of  a 
praanro  variation  an  boat  traaafar  la  mall.  Slqpla  relations  have  baaa  oollactad 
or  daaalapad  ahiah  par  ait  tha  oalaalatim  of  friction  faotora,  recovery  faotore, 
aad  baat  transfer  oaaffioiaata  for  laaiaar  boaadary  layer  flow  aloac  aarfaoaa  of 
oonatant  t— poratnra  with  an  aooaraoy  bat  tor  than  3  par  eant.  A  greater  moor- 
taluty  in  tha  prediction  of  tbo  above  paranatara  ia  oanaad  by  the  poor  knowledge 
of  Fraadtl  aadoer  valnea  for  air  at  high  tanporatnraa.  A  daoldod  advantage  af  the 
ralationa  praaantad  in  tbia  report  lioa  in  the  fact  that  they  appear  to  bold  for 
aay  reasonable  Frandtl  nadOar  variatian. 

halations  far  friction,  raoovary,  and  hast  traaafar  bare  also  barn  reeaanaaded 
far  tatalat  beandary  layers.  A  natbod  bra  bam  daaeribad  by  abiob  tha  off  oat  of 
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oiiCTION  I 


II  .VF»0  i-U  C  TIOi'. 


The  purpose  of  the  present  report  is  to  present  a  survey  and  suiaaary  of 
presently  available  information  on  heat  transfer  at  high  supersonic  speeds.  It  is 
hoped  that  the  paper  is  prepared  in  such  a  way  that  it  is  useful  for  the  designer 
who,  in  order  to  calculate  temperatures  in  aircraft  at  high  speeds,  has  to  have 
information  on  the  convective  heat  transfer  between  the  aircraft  skin  and  the 
atmosphere.  It  is  well  known  that  by  the  aerodynamic  heating  effect  aircraft  which 
fly  with  supersonic  speeds,  are  heated  up  to  such  an  extent  that  the  skin  temperature 
and  the  temperature  inside  the  aircraft  require  the  special  attention  of  the  designer. 

Corresponding  to  the  growing  importance  of  cooling  problems  in  aeronautics,  tne 
number  of  papers  which  deal  with  high  velocity  flow  has  steadily  increaseu  during 
recent  years.  This  is  evidenced  by  the  fact  that  to-uay,  for  example,  each  issue 
of  the  "Journal  of  tne  Aeronautical  sciences*  contains  at  least  one  contribution 
to  the  heat  transfer  problem.  In  this  situation  a  summary  has  to  be  restrictea  in 
its  scope.  In  tnis  report  only  convective  neat  transfer  will  be  discussed  ana 
this  only  for  surfaces  along  wnich  the  pressure  is  constant.  Tne  literature 
collected  in  the  appendix  is  thought  to  be  complete  in  this  field,  ouch  a  restric¬ 
tion  is  permissible  since  aircraft  flying  through  the  atmosphere  at  high  speeus  will 
always  have  very  slender  snapes  so  that  the  pressure  variation  along  the  surface  is 
not  large  and  its  effect  on  heat  transfer  is  of  minor  importance.  Jame  remarks  on 
the  influence  of  pressure  variation  are  however  included,  excluded  also  from  tne 
discussion  will  be  surfaces  which  are  cooled  ty  a  process  known  as  transpiration 
cooling.  In  this  process  the  skin  is  manufactured  from  a  porous  material  and  a 
cooling  medium  is  blown  through  the  skin  into  the  atmosphere  and  builds  up  a  cool 
film  which  protects  the  surface  from  the  influence  of  the  hot  boundary  layer. 
Information  on  this  cooling  process  and  its  influence  on  heat  transfer  under  high 
velocity  conditions  is  still  meager  (Ref.  1,  2,  J ,  4) ■  The  discussion  is  also 
restricted  to  heat  transfer  under  steady  state  conditions.  It  has  been  shown  that 
for  rutes  of  change  in  speed  as  they  are  usually  encountered  In  aircraft  the 
convective  neat  transfer  may  be  regarded  as  quasi  stationary  (See  Appendix  III.)* 

This  means  the  heat  transfer  at  each  instant  may  be  calculated  with  formulae  valid 
for  steady  state  into  which  velocity  and  property  values  are  introduced  as  they  are 
encountered  at  this  specific  instent  (Ref.  3»  6)» 

The  essential  features  of  the  temperature  distribution  as  it  is  encountered 
near  a  wall  moving  with  high  speed  may  be  discussed  with  the  help  of  figure  1.  In 
this  figure  the  temperatures  are  plotted  over  the  distance  from  the  wall.  In  any 
high  speed  gas  flow  two  different  temperatures  have  to  be  distinguished:  the 
static  temperature  T  as  it  would  be  measured  by  an  instrument  which  is  at  rest 
relative  to  the  flow  at  each  location,  and  the  total  Temperature  T.  measured  by  an 
instrument  which  slows  the  flow  adiabatically  down  to  the  veloc ity*zero ,  measured 
relative  to  the  wall.  The  full  lines  represent  the  temperature  distribution  for 
the  case  that  no  heat  is  conducted  from  the  surface  into  the  interior  .. the  wall 
or  transferred  to  the  environment  by  radiation.  The  profiles  of  both  temperatures 
are  shown  in  the  figure. 

TMe  static  temperature  increases  throughout  the  boundary  layer  in  the 
direction  toward  the  wall  because  of  the  heat  generated  by  internal  friction. 

Under  the  conditions  as  specified  no  heat  flow  from  the  wall  into  the  fluid  occurs 
and  therefore  the  temperature  gradient  at  the  wall  is  zero.  However,  tne  wall 
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temperature  is  higher  than  the  static  temperature  in  the  stream  outside  the  boundary 
layer  because  of  the  effect  of  "aerodynamic  heating". 

It  is  interesting  to  observe  the  behavior  of  the  total  temperature  throughout 
the  boundary  layer.  This  temperature  measures  the  total  energy  content  of  the 
fluid,  internal  and  kinetic,  and  this  enerQf  content  aoes  not  change  in  a  steady 
flow  from  which  energy  is  not  extracted,  either  as  heat  or  as  mechanical  energy 
Both  conditions  are  fulfilled  by  the  boundary  layer  as  a  whole.  Therefore,  the 
average  value  of  the  total  temperature  has  to  be  the  same  within  the  bounaary  layer 
as  in  the  outside  flow,  v/hen  a  defect  of  the  total  temperature  in  the  layers  near 
the  wall  is  encountered  for  fluids  with  a  Prandtl  Number  smaller  than  one,  then  it 
has  to  be  compensated  by  excess  of  this  temperature  in  the  outer  region  of  the 
boundary  layer.  For  fluids  with  Prandtl  number  larger  than  one,  the  conditions  ore 
reversed,  the  total  temperature  being  larger  near  the  wall  and  smaller  in  the  outer 
region  of  the  boundary  layer  than  in  free  stream. 

The  difference  between  the  total  and  static  temperatures  at  any  point  is 
determined  by  the  local  velocity  V.  for  a  fluid  with  constant  specific  heat  cp 
this  difference  is  given  by  the  expression  shown  in  the  diagram.  For  variable 
specific  heat  the  difference  between  the  corresponding  enthalpies  i^  •  i  is  equal 
to  the  kinetic  energy  of  the  air  V*/2.  For  gases  under  the  condition  of  no  heat 
transfer,  this  wall  temperature  which  will  be  called  in  this  report  the  "recovery 
temperature"  is  almost  equal  to  the  total  temperature  in  the  stream  outside  the 
boundary  layer.  This  leads  at  supersonic  velocities  to  very  high  wall  temperatures, 
fortunately,  conditions  are  usually  much  more  favorable  on  the  skin  of  aircraft 
flying  at  high  speed  because  heat  is  radiated  a way  from  the  surface  of  the  aircraft 
or,  for  short  flight  duration,  heat  is  being  stored  in  the  skin  or  other  structural 
elements.  The  temperature  profile  found  in  the  boundary  layer  under  such  circum¬ 
stances  is  shown  by  dashed  lines  in  figure  1.  It  is  interesting  to  note  the 
peculiar  shapes  which  the  total  temperature  profile  may  assume.  This  profile  is 
the  one  which  can  be  comparatively  easily  measured  by  total  temperature  probes, 
whereas  today  no  satisfactory  direct  measurement  of  local  static  temperatures  is 
possible. 

from  the  discussion  of  figure  1  it  is  apparent  that  heat  transfer  in  high 
speed  flow  is  essentially  influenced  by  two  parameters  which  are  usually  not 
present  in  other  heat  transfer  problems i  a  marked  effect  of  frictional  heating 
and  a  large  temperature  variation  which  occurrs  in  the  neighborhood  of  surfaces  in 
a  high  velocity  flow  field,  even  when  the  surface  itself  is  cooled  down  to  the 
outside  stream  temperature. 

It  will  be  shown  later,  that  for  the  case  when  the  wall -temperature  is  cooled 
to  the  stream  static  temperature  a  maximum  temperature  within  the  boundary  layer 
arises  which  for  a  I-ach  number  10  is  five  time  the  stream  temperature  or,  for  a 
i-ach  number  20  is  equal  to  18  tij.es  the  stream  temperature.  These  temperature 
variations  are  unusually  large  and  cause  the  property  values  to  change  extensively 
throughout  the  boundary  layer.  The  temperature  within  the  boundary  layer  may 
become  so  large  that  it  causes  dissociation  of  the  air. 

frictional  heating  as  determined  by  the  velocity  of  the  air  stream,  large 
temperature  variation,  and  high  temperatures  are  the  parameters  which  will  have  to 
be  discussed  extensively  in  this  report  with  regard  to  their  influence  on  heat 
transfer , 
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mother  condition  which  is  connected  v;ith  aircraft  application  is  tne  lou 
. ressure  ana  correspondingly  tie  low  density  of  toe  air,  encountered  by  tue  eircreft 
at  hijh  flight  altitude.  Tne  air  at  ni^pi  ultitu-es  is  so  rarified  taut  the  lengths 
of  the  paths  'which  t.io  single  air  molecules  travel  between  collisions  became 
considerable  and  of  the  suae  order  of  ua^yiitu^e  as  the  thickness  of  the  boundary 
layer .  Ibis  cn-n^es  heat  transfer  conditions  from  what  is  usually  experienced  in 
air  flow  under  normal  pressure.  This  factor  will  ulso  be  discussed  in  its 
influence  on  neat  transfer  in  one  of  the  following  chapters. 


oJCTICI.  II 
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Jectijn  III  will  deal  extensively  witii  tue  frioticn  and  neat  transfer 
cant.itijn3  aa  t  ie;-  exist  ij  lcuincr  or  turbulent  boundary  layers  on  surfaces  in 
two-i L.ensional  flow  for  which  t.ie  pressure  does  not  change  in  flow  direction. 

Juc.i  conditions  arc  realized  on  a  flat  plate  parallel  to  the  stream  direotion  as 
long  as  the  boundary  layers  are  very  thin  so  that  the  displacement  of  the  flow  by 
this  boundary  layer  can  be  neglected. 

Constant  pressure  is  in  supersonic  flow  also  encountered  along  the  surface  of 
wedges  and  cones  when  the  shock  wares  by  which  the  oncoming  flow  is  turned  into  e 
direction  parallel  to  the  surface  are  attached  to  the  vertex  of  the  objeot.  /or 
this  case  simple  conversion  formulas  exist  by  which  friction  and  heat  transfer  data 
for  cones  and  wedges  can  be  obtained  from  information  available  for  flat  plates. 
These  will  be  discussed  in  this  chapter.  duch  conversion  rules  are  important 
because  parts  of  missiles  or  aircraft  have  conical  shapes  end  also  because  many 
test  results  are  available  on  cones  and  wedges  and  can  in  this  way  be  compered  with 
fist  plate  data. 

friction  and  heat  transi’er  on  the  surface  of  a  cylinder  in  a  flow  parallel  to 
its  axis  are  identical  with  the  ones  on  a  flat  plate  as  long  ae  long  as  the  smallest 
radius  of  the  curvature  of  the  cylinder  cross  section  is  large  compared  with  the 
boundary  layer  thickness,  dome  investigations  as  to  the  limiting  conditions  for 
this  analogy  will  also  be  discussed. 

Far  a  wedge  in  a  supersonic  flow  which  is  directed  perpendicular  to  its  edga 
and  which  has  attached  shock  waves  the  flow  conditions  are  constant  in  the  field 
between  the  shocks  and  the  surface  of  the  wedge.  Constant  i-aoh  number,  constant 
total  and  static  temperatures  and  constant  density  prevail  in  this  field,  dines 
the  flow  on  such  a  wedge  is  two-dimensional  and  since  the  pressure  along  the 
surface  of  the  wedge  is  constant,  the  relationship  for  skin  friction  and  heat 
transfer,  valid  for  flat  plates,  may  be  immediately  applied  to  the  surfaces  of  such 
a  wedge.  Conditions  between  the  shock  and  the  surface  correspond  to  the  flww 
conditions  at  the  flat  plate  outside  the  boundary  layer. 

A  cone  may  now  be  considered  which  is  located  in  a  supersonic  flow  directed 
parallel  to  the  cone  axis.  The  onstream  conditions  may  be  such  that  the  shock  were 
is  attached  to  the  apex  of  the  cone.  The  shock  has  in  this  case  a  conical  ahqpe 
with  axis  and  apex  coincident  with  the  solid  cone.  Stream  conditions  in  this  ease 
ere  not  constant  in  the  field  between  this  shock  end  the  cone  surface,  however, 
they  are  constant  along  any  conical  surface  with  the  same  axis  and  apex  am  long  as 
the  displacement  of  the  flow  by  the  boundary  layer  is  oeglegible.  Jbr  thin  boundary 
layers,  therefore,  this  condition  will  also  hold  along  a  conical  envelope  very  near 
to  the  cone  surface  and  just  outside  the  boundary  layer,  tfhen  friction  or  heat 
transfer  data  from  flat  plates  ere  to  be  converted  to  that  on  the  surfaces  of  eones 
then  the  flow  conditions,  existing  in  this  envelope  which  is  praetieclly  coincident 
with  the  cone  surface,  correspond  to  the  flow  conditions  for  the  plate,  ^ain  the 
pressure  is  constant  along  the  surface  in  flow  direction. 
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The  flow,  however,  is  not  two-dimensional  anti  therefore  conversion  formulas 
have  to  be  used  in  order  to  obtain  from  flat  plate  data  relationships  valid  for  the 
surface  of  cones.  Hantzsche  and  Wendt  (Ref.  12)  show  by  a  transformation  of  the 
laminar  boundary  layer  equations  for  flow  along  the  cone  surface  to  the  flat  plate 
case  that  the  temperature  recovery  factor  remains  unchanged  fay  this  transformation 
and  has  the  same  value  for  cones  as  for  flat  plate.  Local  friction  coefficients 
and  heat  transfer  coefficients,  valid  for  the  flat  plate,  have  to  be  multiplied  fay 
V3T  to  obtain  friction  and  heat  transfer  coefficients  on  the  surface  of  the  cons  at 
the  same  distance  from  the  leading  edge  and  under  identical  stream  conditions. 

This  fact  may  also  be  stated  in  a  different  way,  namely,  that  local  friction  and 
heat  transfer  coefficients  are  id  ntical  for  cones  and  flat  plates  ut  locations  for 
which  the  cone  Reynolds  number  is  times  the  value  of  the  flat  plate  Reynolds 
number.  At  such  corresponding  locations  within  the  boundary  leyer  arc  completely 
identical  and  therefore  not  only  the  heat  transfer  coefficient  and  friction 
factors  but  also  the  stubility  conditions  within  the  boundary  layer  and  therefore 
the  conditions  for  transition  to  turbulence  are  the  same. 

An  analogous  transformation  on  the  turbulent  boundary  layer  equations  written 
in  time  mean  values  has  been  made  by  Gazley  (Ref.  11),  and  in  a  veiy  general  font 
taking  into  account  the  variation  of  the  property  values  by  Van  Driest  (Ref.  9). 

The  result  of  this  transformation  is  that  again  conditions  within  the  boundary 
layer  are  identical  on  such  locations  at  the  flat  plate  and  on  the  cone  surfaoe  for 
which  the  cone  Reynolds  number  is  equal  to  twice  the  flat  plate  Reynolds  number. 

%  the  sikple  rules  stated,  friction  and  heat  transfer  can  be  calculated  from 
formulas  which  are  valid  for  flat  plates. 

For  some  purposes  It  is  important  to  know  also  how  much  the  local  and  the 
average  friction  and  heat  transfer  values  differ  on  the  cone  surface  and  the  flat 
plate  at  locations  with  the  same  Reynolds  number.  For  laminar  flow  it  is  already 
stated  that  the  local  heat  transfer  coefficient  and  friction  factor  on  the  cone  are 
by  the  factor  VT  larger  than  on  the  flat  plate.  The  average  heat  transfer  coeffi¬ 
cient  and  friction  factor  for  the  cone  are  larger  by  the  factor  2/  ^  4  V.  It 
may  be  seen  that  the  average  friction  and  heat  transfer  values  do  not  differ  much 
for  canes  and  flat  plates  when  they  are  compared  at  the  same  Reynolds  number.  For 
the  turbulent  friction  coefficient,  using  Blsslus  equation  and  Reynolds  analogy, 
the  following  relationships  for  a  turbulent  boundary  layer  can  be  derived:  The 
local  friction  and  heat  transfer  coefficients  are  by  the  factor  2***  1.15  larger 
than  for  a  flat  plate.  For  th->  average  values  the  factor  is  1.022.  Therefore, 
the  difference  between  average  friction  and  he*t  transfer  coefficients  on  cones 
and  flat  plates  considering  locations  with  the  same  Reynolds  number  is  only  2  per 
cent  for  turbulent  flow.  It  may  be  added  here  that  a  transformation  of  the  boundary 
layer  equations  for  arbitrary  bodies  of  revolution  in  rotationally  synrcetric  flow 
to  the  two-dimensional  case  has  be*>n  presented  in  Reference  14* 

For  a  cylinder  with  circular  cross  section  in  a  flow-subsonic  or  ougsvaeaio— 
which  la  parallel  to  the  cylinder  axis  conditions  within  the  boundary  layer  are 
principally  different  froi.i  the  conditions  on  a  flat  plate  uinoe  cylindrical  areas 
within  the  boundary  layer  parallel  to  the  surface  increase  in  area  proportional  to 
the  distance  from  the  cylinder  uxi3.  flic  influence  of  this  condition  on  heat 
transfer  end  friction  in  a  turbulent  boundary  layer  has  been  determined  analytically 
fay  Jakob  and  Low  (Ref.  13),  by  it.  U.  Lckert  (Ref.  10),  und  by  Seban  (Ref.  17).  lbs 
latter  investigator  found  that  thin  condition  does  not  influence  the  recovery  factor 
so  that  flat  plate  recovery  factors  can  be  iumediately  applied  to  cylinders  in  axial 
flow. 
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Jith  regard  to  the  average  skin  friction  the  different  predictions  vary 
widely.  An  experimental  investigation  in  reference  65,  mude  on  cylinder-cone 
combinations,  found  that  at  a  Reynolds  number  around  10  an  increase  of  the  average 
friction  factor  of  5  percent  was  connected  with  a  change  of  the  length  L  to  diameter 
d  ratio  from  0  to  23.  The  calculation  by  Jakob  and  Dow  gives  for  a  change  of  l/d 
from  8  to  23  an  increase  of  15  percent.  H.  U.  ’’fekert’s  method  results  for  the  same 
conditions  In  1.5  percent  increase.  Th»  evidence  presented  does  not  give  a  conclu¬ 
sive  quantitative  information  on  how  much  friction  and  heat  transfer  on  a  cylinder 
are  different  from  the  ones  on  a  flat  plate.  It  can  be  concluded  hov“v*r,  that  up 
to  a  ratio  of  boundary  layer  displacement  thickness  to  cylinder  radius  of  1:100  to 
1:50,  the  differenoe  between  cylinder  and  flat  plate  data  should  be  negligible. 

For  cylinders  in  a  flow  normal  to  their  axis  the  condition  for  heat  transfer 
ere  rather  involved,  a  boundary  layer  builds  up  around  the  upstream  part  of  the 
surface.  The  large  pressure  variation  along  the  surface  makes  the  development  of 
this  boUhdary  layer  Quite  different  from  the  conditions  considered  before  in  this 
chapter.  Hie  effect  of  a  pressure  variation  on  the  boundary  layer  will  be  considered 
in  somewhat  more  detail  in  a  later  chapter.  Heat  transfer  to  the  downstream  part  of 
the  cylinder  where  the  flow  separates  from  the  surface  is  still  poorly  understood. 
Extremely  low  recovery  temperatures  (partially  below  static  upstream  -  temperature) 
have  been  observed  in  this  region  at  subsonic  upstream  speeds.  (Ref.  176.  179) 
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aaCTXON  III 


HLftT  TRaKo/JlR  TO  k  FUT  FLkTh  ,/ITH  CONSTAT  SURFACE  TU>JPUUTURB 


This  section  deals  with  what  ia  usually  called  the  flat  plate  case.  More 
exactly,  it  deals  with  heat  transfer  in  two-dimensional  flow,  the  pressure  of  which 
is  constant  along  the  surface,  for  thin  boundary  layers  this  type  of  flow  is  found 
on  e  flat  plate.  There  are,  however,  special  conditions  like  very  low  pressure, 
very  high  velocity  where  boundary  layers  become  quite  thick.  The  flow  outside  the 
boundary  layer  is  then  displaced  and  curved  along  a  flat  plate,  according  to  the 
increase  of  the  displacement  thickness  of  the  boundary  layer  in  downstream  direction. 
As  a  consequence  of  this,  the  pressure  decreases  in  flow  direction.  Constant 
pressure  in  such  a  case  would  be  rather  obtained  by  a  slightly  concave  surface.  In 
supersonic  flow,  constant  pressure  along  a  surface  is  ulso  found  on  welges  and  conns 
with  attached  shocks.  The  simple  rules  with  which  flat  plate  formulas  for  heat 
transfer  and  friction  can  be  transferred  to  the  conditions  on  wedges  and  cones  have 
been  discussed  in  the  previous  chapter. 

a  constant  temperature  along  the  surface  is  also  specified  for  this  ohapter . 
Generally,  on  missiles  and  aircraft  the  temperature  will  vary  along  the  surfaoe 
and  it  is  known  today  that  such  a  temperature  variation  has  a  considerable  effect 
on  heat  transfer.  The  plate  with  constant  temperature,  however,  serves  as  a 
standard  condition  to  which  comparatively  simple  formulas  apply.  A  variation  of 
the  surfaces  temperature  has  to  be  dealt  with  specifically  by  calculations  in  each 
individual  case.  This  procedure  will  be  discussed  in  a  later  ohapter. 

Conditions  on  the  wall  surface  of  the  plate  will  be  denoted  by  a  subscript  w, 
conditions  in  the  stream  outside  the  boundary  layer  by  a  subscript  s.  Harm  it  is 
assumed  that  flow  conditions  in  the  field  outside  the  boundary  layer  (velocity  Yg. 
temperature  Tg,  density  #8)  are  constant.  The  symbol  T  denotes  absolute  taajeratarea 
(in  °H. ) .  In  the  flow  T  means  always  static  temperature.  The  total  temperature 
(as  measured  by  a  total  temperature  probe  at  rest  relative  to  the  plate)  is 
indicated  by  a  subscript  t.  The  following  dimensionless  flow  parameters  are  useds 


(B  -  1) 


1, 


(B  -  2) 


The  Reynolds  number  Re  is  based  on  the  distance  X  from  the  leading  edge  of  the 
plate .  The  property  values  p  and  Jn  in  the  Reynolds  number  will  be  introduced  at 
different  locations  (wall  or  stream),  whereas  the  hash  number  Ms  is  always  based  on 
the  sound  velocity  as  at  stream  temperature. 
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The  local  shearing  stress  w  which  the  flew  exerts  on  the  plate  surface  is 
expressed  by  a  dimensionless  friction  factor  defined  in  the  following  way 

_g 

ft-t  (B  -  3) 

The  local  convective  heat  flow  per  unit  time  and  area  at  the  plate  surface  is 
described  by  a  heat  transfer  coefficient  h  which  is  defined  in  the  following  way 

%  -  h  (fr  -  \  )  (B  -  4) 


T  is  the  actual  wall  temperature  and  Tr  is  the  temperature  which  the  speoifie 
location  on  the  surface  assumss  when  the  oonvective  heat  transfer  is  zero.  This 
temperature  will  be  called  ’recovery  temperature*  since  it  is  conventionally 
described  by  the  temperature  recovery  factor  Sq.  (Ba-4).  The  heat  transfer 
coefficient  is  in  turn  expressed  non-dime naionally  by  the  Nusselt  number 


In  - 


(B  -  5) 


or  by  the  Stanton  number 


(B  -  6) 


UKURbR  BOUNDARY  LAY'S  FLOW 


The  differential  equations  which  describe  the  flow  and  energy  conditions  for  a 
laminar  boundary  layer  are  well  established  under  normal  pressure  conditions  la  the 
continuum  flow  regime. »  Numerous  solutions  of  these  differential  equations  have  beam 
worked  out  and  are  today  available  to  give  information  in  a  wide  range  of  conditions, 
extending  for  instance  to  very  high  Mach  number  and  large  temperature  differences. 
These  results  generally  check  satisfactorily  with  experimental  investigations,  whan 
the  experiments  were  made  at  the  boundary  conditions  which  correspond  to  the  boundary 
layer  eolations.  Therefore ,  in  the  laminar  range,  the  most  reliable  information  sen 
be  obtained  from  calculations .  The  main  difficulty  in  solving  the  differential 
equations  is  connected  with  the  fact  that  the  property  values  (viscosity,  heat 
conductivity,  density,  and  specific  heat)  vary  for  all  gases  mad  fluids  with 
temperature  and  pressure.  It  is  quite  difficult  to  take  into  proper  account  the 
actual  variation  of  the  property  values  in  the  solution  of  the  differential 
equations.  The  flat  plate  case  is  in  this  reapeet  simpler  than  flow  end  heat 
transfer  for  general  conditions.  At  the  flat  plats,  as  mentioned  above,  the 
pressure  is  constant  along  the  surface.  In  the  thin  boundary  layer  the  pressure 
also  does  not  change  in  e  direction  perpendicular  to  the  surface  so  that  it  is  then 
constant  over  the  whole  field  in  which  the  flow  end  energy  transfer  occur.  In 
consequence  of  this  fact  only  the  variation  of  density  with  temperature  has  to  be 
accounted  for  in  solutions  of  tbs  boundary  layer  equations  for  the  flat  pints.  In 
this  connection  it  night  be  pointed  out  that  the  terminology  which  has  besoms 
customary,  namely  to  cell  a  boundary  layer  on  a  flat  plate  which  develops#  under 
high  flow  velocities  e  'compressible  boundary  layer*  is  not  too  well  chosen,  the 
established  meaning  of  the  term  in  physics  is  e  change  of  volume  in  oonaequenos  of  a 
Chengs  in  pressure  end  such  e  change  does  not  occur  in  a  boundary  layer  on  a  flat 
plate  since  the  pressure  in  the  whole  field  of  interest  is  constant.  This  expresses 
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itMlf  in  the  fact  that  the  character  of  the  boundary  layer  equation*  ia  exactly 
the  same  for  subsonic  as  in  the  supersonic  range. 

Constant  Property  Solutions 

The  first  solutions  of  the  laminar  boundary  layer  equations  were  obtained  with 
the  assumption  that  the  property  values  in  these  equations  are  constant  (independent 
of  pressure  and  tec.perature ) .  Under  this  condition  the  local  skin  friction  factor 
was  found  to  be 


cf  -  0.6M /Vti 


(Be  -  1) 


The  local  Nusselts  number  was  determined  as 

la  -  0.332  >fe  far 


(Ba  -  2) 


when  the  heat  transfer  coefficient  h  on  which  the  NUsselt  number  is  baaed  is 
defined  by  the  following  equation 


S  -  Mr,-*.) 

(B*  -3) 

in  which  Tr  ia  the  temperature  which  the  plate  assumes  when  it  is  only  in  convective 
heat  exchange  with  the  high  speed  flow  but  otherwise  not  cooled  or  heated  (by 
radiation,  or  by  heat  conduction  into  the  solid  wall  interior).  This  wall  teeqpera- 
ture  is  ia  high  speed  flow  larger  than  the  static  temperature  in  the  stream  outside 
the  boundary  layer  because  of  the  aerodynamic  heating  effect  within  the  boundary 
layer  and  its  value  has  to  be  known  if  the  heat  flow  on  a  surface  under  high  spaed 
conditions  is  to  be  calculated.  The  calculations  which  wars  made  under  the 
assumptions  of  constant  property  values  showed  that  the  following  temperature  ratios 


lull 

which  is  called  temperature  recovery  factor  is  a  function  of  only  the  Drandtl 
number  of  the  fluid.  This  function  is  shown  over  a  wide  range  of  Frandtl  numbers 
in  Tleur*  2.  (Ref.  37) •  It  can  ba  approximated  by: 


r  -  yff 


(B»  -  5) 


in  the  range  of  Frandtl  numbers  between  0.5  end  5.  The  wall  temperature  T  which 
is  deaerlbed  by  the  temperature  recovery  factor  will  in  this  report  be  called 
"recovery  temperature" . 


9 


The  heat  transfer  coefficient  is  often  expressed  in  a  dimensionless  form  by  s 
parameter  celled  otanton  number: 


St  «  h  /  f-Cp-T,  (fia  .  6) 

Zt  can  easily  be  seen  that  the  following  connection  exists  between  the  Nusselt 
number  and  the  Stanton  Number: 


•  Mo  /  Ie*ft» 


(Ba  -  7) 


A  comparison  of  Equation  Ba  -  1  for  the  skin  friction  factor  and  Equation  Ba  -  2 
for  the  local  ftuseelt  number  results  in  the  following: 


St  -  0.5  ef  Pr"**  (Ba  -  8) 

This  relation  between  skin  friction  and  heat  transfer  is  useful  because  it  still 
holds  under  certain  conditions  when  the  property  ▼  slues  vary  with  temperature.  This 
will  be  discussed  in  the  next  paragraph. 

Reynold's  analog  between  skin  friction  and  heat  transfer  can  be  expressed  by 
the  following  equation: 


The  heat  transfer  calculated  from  Reynold's  analogy  therefore  differs  from  the 
correct  relationship  for  laminar  flow  by  the  factor  For  gases  the  FTandtl 

number  is  not  too  different  from  1  and  therefore  Reynolds  analogy  gives  a  reasonable 
first  approximation  to  the  real  heat  transfer  value.  The  relationships  which  have 
been  detersdned  for  constant  property  values  can  be  expected  to  hold  satisfactorily 
for  real  fluids  or  gases  as  long  as  the  temperature  variation  throughout  the 
boundary  layer  is  comparatively  small.  It  has  bean  found  that  this  is  the  coos  as 
long  as  the  ratio  of  wall  to  static  fires  stream  temperature  is  not  too  large  (say 
between  0.5  and  1.3)*  Property  data  for  determining  the  Reynolds  and  Brudtl 
austere  may  be  formed  with  temperatures  or  T_.  Hie  advantage  of  the  constant 
property  solutions  is  that  they  depend  on  e  minimum  of  parameters  and  that  they 
apply  generally  to  any  fluid. 

fc^itlona  For 

In  the  eupersonio  range  and  for  large  temperature  differences  tbs  feet  that 
the  property  values  actually  depend  on  temperature  has  to  be  accounted  for  in  a 
solution  of  the  laminar  boundary  lqrer  aquations,  figure  3  shows  that  the  viscosity 
and  heat  conductivity  vary  considerably  faster  with  temperature  than  the  specific 
heat  and  the  Frond  tl  nurber.  Consequently  the  next  step  wee  to  obtain  boundary 
layer  solutions  for  variable  viscosity  end  conductivity  but  constant  Rrandtl  r 
end  spec if ic  beet.  The  equation 


t 


Zt 

•p 
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(Be  -  10) 


si;3..s  tl:-t  -'or  constant  lr^’ic'tl  number  end  apecii'ic  heat,  the  heat  conductivity 
varies  proportion.,  liy  to  the  TiModtf.  Therefore  the  law  of  temperature  aepeneence 
lies  to  be  prcacritei.  for  one  property  only.  The  following  law  was  assumed  to 
e.cpress  the  te:  .peraturc  dependency  of  the.  vicosity 


A 

A> 


(a»  -  ii) 


ILand  ?Q  ere  the  viscosity  and  temperature  at  some  reference  condition.  The  uf-  of 
this  law  brings  two  advantages,  first  it  simplifies  the  calculation  procedure  a.  a 
secondly  it  melees  the  result  dependent  on  a  characteristic  temperature  ratio  only, 
for  instance,  on  the  ratio  of  wall  temperature,  to  free  stream  static  or  total 
temperature .  This  temperature  ratio  as  well  as  the  exponent  H  have  to  be  added  to 
the  parameters  mentioned  above  for  the  constant  property  value  solutions. 


Equation  Be  -  11  expresses  the  aetual  temperature  variation  of  the  vicosity 
for  gases  t;nd  specifically  for  air  reasonably  well  in  a  temperature  range  which  is 
not  too  large.  However,  it  is  found  that  the  numerical  value  of  the  exponent 
depends  on  the  temperature  level,  for  very  low  temperatures  a) is  approximately  1 
and  it  decreases  with  increasing  temperature  asymptotically  towards  the  value  0.5. 
Consequently,  in  different  calculations  different  numerical  values  for  this 
parameter  were  used.  The  following  general  results  were  obtained:  (for  a  summary 
of  these  calculations  see  Reference  45) •  It  was  found  that  the  expression  Be  -  5 
for  the  recovery  factor  still  holds  under  the  assumptions  of  this  paragraph.  In 
the  same  way  the  relationship  Ba  -  8  between  akin  friction  and  heat  transfer  is 
still  valid.  On  the  other  hand  it  was  noted  that  the  akin  friction  and  the  hast 
transfer  vary  now  with  taeh  number,  with  the  characteristic  temperature  ratio,  with 
the  value  of  the  parameter  41  ,  and  the  temperature  level.  Only  the  variation  of  the 
friction  coefficient  has  to  be  discussed  because  of  the  validity  of  equation  Be  -  8. 
Specifically  the  variation  of  the  skin  friction  coefficient  with  the  mentioned 
parameters  depends  on  the  way  in  which  they  are  defined.  Aram  equations  B  -  1  and 
E  -  3  which  define  the  skin  friction  ooeffioient  end  the  Reynolds  number  it  can  be 
seen  that  for  variable  property  values  a  specific  temperature  must  be  chosen  at 
which  these  values  are  introduced  into  the  right  hand  side  of  the  aquations. 

Usually  the  free  stream  static  temperature  ia  chosen  as  such  a  temperature .  In 
this  case  the  akin  friction  coefficient  is  found  to  decrease  considerably  with 
increasing  i-nch  number,  with  increasing  temperature  ratio  Ty/Ts,  and  with  decreasing 
value  of  the  parameter  Ii  .  figure  4  taken  from  reference  143  shows  the  dependence 
of  the  friction  factor  on  Reoh  number.  The  dashed  curves  and  the  calculations  by 
Van  Ikies t,  and  Young  and  Janssen  will  be  discussed  later  on.  It  is  Interesting  to 
note  that  a  specific  value  4)  ■  1  makes  the  dependence  of  friction  factor  on  fern  and 
T  /t  vanish.  The  relations  Ba  -  1  to  Ba  -  9  add,  therefore,  not  only  for  constant 
properties  but  also  when  viscosity  and  hast  conductivity  vary  proportionally  to  the 
absolute  temperature.  This  fact  adds  considerably  to  the  importance  of  the  constant 
property  solutions. 
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If  the  property  values  in  Equations  B  -  1  and  B  -  3  for  the  skin  friction 
factor  for  *1  different  from  one  are  introduced  at  the  wall  temperature,  then  the 
trend  of  the  variation  of  the  skin  friction  factor  with  hach  number  is  opposite  to 
the  one  mentioned  above. 

It  may  therefore  be  expected  that  property  values  introduced  at  a  proper 
reference  temperature  situated  somewhere  between  the  temperature  extremes  encountered 
within  .he  boundary  layer  will  cause  the  variation  of  the  friction  factor  with  to 
and  Ty/Ta  to  vanish,  actually,  Rube sin  and  Johnson  (Ref.  49)  have  shown  that  it  is 
possibleto  find  such  a  reference  temperature  T*  and  that  it  can  be  expressed  with 
good  accuracy  by  the  following  equation  taken  from  Rubesin's  and  Johnson's  paper  1 

r  -  1  ♦  O.032  (Ik)2  ♦  0.58  (fr  -  1) 

•  •  (Ba  -  12) 

This  expression  for  a  reference  temperature  can  also  be  transformed  into  the 
following: 

t*  -  T,  -  0.58  (%  -  T.)  ♦  0.19  (Tr  -  T,) 

(Be  -  13) 

In  this  form  the  relation  for  the  reference  temperature  appears  especially 
instructive,  as  can  be  seen  for  same  special  oases  sketched  in  figure  3: 

ftor  low  flow  velocities  the  recovery  wall  temperature  T_  is  practically  equal 
to  the  stream  temperature  Tg.  The  temperature  profile  in  this  case  has  the  shape 
indicated  in  the  upper  diagram  of  the^figure  and  Equation  Ba  -  13  shows  thut  in 
this  case  the  reference  temperature  T  is  located  approximately  hulfway  between 
wall  temperature  and  stream  temperature  Tg.  This  agrees  with  the  results  of 
other  investigations  on  heat  transfer  at  low  velocities. 

to r  high  velocities  but  a  small  rate  of  heat  transfer  the  wall  temperature 
Tw  is  approximately  equal  to  the  reoovery  temperature  Tr.  The  temperature  profile 
is  for  this  case  indicated  in  the  center  portion  of  the  figure  and  Equation  Ba-  13 
shows  that  the  difference  between  reference  temperature  T*  and  stream  temperature 
Ts  le  now  77%  of  the  difference  between  wall  temperature  %  and  stream  temperature 
Tg.  Jince  the  fluid  layers  near  the  wall  will  Influence  the  heat  transfer  to  the 
wall  moat  strongly  and  since  these  layers  have  now,  in  an  average,  a  higher 
temperature  than  in  the  first  case  considered,  the  value  for  the  reference 
temperature  is  very  credible. 

Am  a  third  example  let  us  assume  that  the  wall  is  eooled  to  the  stream 
temperature  as  indicated  in  the  lower  diagram  of  the  figure.  It  will  be  seen 
later  on  that  in  this  case  the  difference  between  the  maximum  temperature  T_ 
within  the  boundary  layer  and  the  stream  temperature  is  approximately  25%  of  the 
difference  between  recovery  temperature  and  stream  temperature,  tor  the  difference 
between  reference  temperature  and  stream  temperature  the  corresponding  value  la 
19X,  which  is  again  possible.  It  should  be  pointed  out  that  figure  4  presents 
friction  factors  for  the  case  where  the  wall  is  at  recovery  temperature,  fbr  a 
cooled  wall  the  friction  factors  ere  higher  and  for  a  heated  wall  lower  than  the 
values  in  the  figure.  Calculations  of  the  friction  faotor  with  the  reference 
temperature  as  given  by  Equations  Ba  -  12  and  Be  -  13  properly  accounts  for  these 
variations  es  will  be  shown  in  the  following  chapter. 
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Solutions  for  Air  with  Actual  Temperature  Variation  of  All  Properties 


For  very  high  velocities  as  are  expected  in  future  aeronautical  applications 
the  temperature  variation  thoroughout  the  boundary  layer  becomes  so  large,  even  in 
the  cese  where  the  well  la  cooled  down  to  practioally  the  free-atreem  static 
temperature ,  that  it  is  doubtful  how  valid  the  results  of  calculations  are  which  are 
based  on  the  assumptions  mentioned  in  the  proceeding  paragraph.  Accordingly,  a 
number  of  recently  published  papers  oonsider  the  specific  heat  and  the  frandtl 
nun.be r  as  well  as  the  heat  conductivity  and  the  viscosity  to  vary  with  temperature . 
all  of  the  calculations  have  been  made  for  air  and  more  accurate  relationships 
than  a  power  law  for  the  variation  of  the  vicosity  and  the  heat  conductivity  with 
temperature  have  been  ussumed. 

Sutherland's  relationship: 

%  ■ 

(flu  -  14) 

with  the  constant  C  *  392°  R.  for  air,  or  tabulated  values  are  used  to  express  the 
viscosity*  Tabulated  values  are  used  for  the  specific  heat  and  lxandtl  number. 

The  conductivity  is  again  fixed  by  equation  Bu  -  10.  The  main  difference  in  the 
various  sets  of  calculations  is  caused  by  the  fact  that  differeht  laws  for  the 
Brandtl  number  dependence  on  temperature  have  been  used.  Klunker  and  lacLean 
(Ref.  41)  us  well  as  Young  and  Janssen  (Aef .  J>9)  use  Prandtl  number  data  as  contained 
in  the  Gas  Tables  and  in  a  paper  by  T*ibua  and  fioelter  (Ref.  177)*  Van  driest 
(Ref.  36)  on  the  other  hand  bases  his  calculations  on  Rrandtl  number  values  which 
have  been  reconnended  by  the  National  haem,  of  .Standard  (Ref.  1£6).  Figure  3 
indicates  the  Rrandtl  numbers  as  used  in  both  sets  of  data.  It  may  be  recognized 
how  large  the  differences  in  both  sets  of  frandtl  number  data  are.  All  of  the 
calculations  hod  to  extrupolute  Rrandtl  number  values  to  larger  temperatures. 

This  is  indicated  by  the  dashed  part  of  the  lines.  The  large  differences  in  both 
sets  of  data  reflect  the  lack  of  aocurate  knowledge  of  heat  conductivity  and 
Rrandtl  number  for  air  at  hitter  temperatures  and  is  caused  by  the  fact  that  it  is 
ext rerely  difficult  to  measure  heat  conductivity  values  for  gases  at  larger 
temperatures ■ 

The  calculations  mentioned  above  reflect  the  real  condition^  as  encountered  fay 
an  aircraft  in  che  atmosphere  better  than  the  ones  discussed  in  the  previous 
paragraph  since  they  are  based  on  more  accurate  property  value  data;  however,  they 
have  one  disadvantage,  namely,  that  flow  ana  heat  transfer  now  appear  to  depend  on 
one  additional  parameter  namely,  on  the  temperature  level  in  the  air  stream.  The 
governing  parameters  are  now  Reynold's  number,  1-mch  number,  Irandtl  number,  a 
characteristic  temperature  ratio  and  a  characteristic  temperature,  for  instance  wall 
to  stream  temperature  ratio,  and  stream  temperature.  The  number  of  parameters  is 
considerable  and  correspondingly  the  calculations  become  rather  extended. 
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All  the  calculation*  show  that  the  frietion  factor  varies  proportional  to  the 
square  root  of  the  Reynold's  number  in  the  eame  way  as  in  Iquation  Be  -  1.  Instead 
of  the  constant  0.664  in  this  equation,  however,  a  function  of  Mach  number,  tempera¬ 
ture  ratio  and  temperature  level  appears.  The  most  extensive  calculations  of  this 
nature  have  been  done  by  Young  and  Janssen  (Ref.  59)*  The  authors  calculated 
frietion  factors  fora  wide  range  of  Mach  numbers  (up  to  3°)  for  a  large  range  of 
wall  to  stream  temperature  ratios  (1.3  to  22)  and  for  three  values  of  the  static 
stream  temperature  (100,400,800°  R.).  They  also  determined  a  reference  temperature 
such  that  introduction  of  the  property  values  at  that  temperature  makes  the  value 
CU  equal  to  the  incompressible  value  O.664.  Young  and  Janssen  reoocsiended 
^nation  Ba  -  12  for  Ms  up  to  5  and  give  the  following  relationships  which 
approximate  the  calculated  reference  temperatures  well  in  a  Mach  number  range  from 
3  to  10> 


fori  %  •  Tr  1  T*A#  -  0.0956  [l 3(V%)  ♦  6]  (Ba  -  13) 
fori  Tv  -  Tri  T*Ag  -  o.TO  ♦  0.023  («a)2  ♦  o.5*  (%A,) 


Since  it  is  inconvenient  to  use  different  equations  for  different  tech  number 
ranges  and  far  the  wall  at  recovery  temperature  and  at  a  different  value,  it  was 
investigated  whether  the  reference  temperatures  can  be  expressed  in  the  whole  Mach 
number  range  by  one  relationship.  The  calculations  presented  in  .Appendix  I 
resulted  in  the  following  equations  for  the  reference  temperature 

1*  •  %  ♦  o.5  <TW  -  f  #)  ♦  o.22  (»r  -  T.) 

(Ba  -  16) 

The  equation  may  also  be  written  as  follows  1 


•  0.5  (T,  ♦  Tw)  ♦  o.22  rt^L  Ufc)2  f, 

(Ba  -  16a) 

with  r  indicating  the  recovery  factor  and  f  the  ratio  of  specific  beets. 

Table  la  contains  in  the  fourth  oolumn  the  friction  parameters  Gg  taken 
from  Ref  39*  In  the  reference  the  property  valuta  were  baaed  on  the  stream 
temperature .  The  fifth  column  oontaina  the  amas  frietion  parameter  as  calculated 
with  the  reference  temperature  given  in  the  last  aquation  in  the  following  way  1 


Aquation  Ba  -  1  was  assumed  to  hold  whan  the  properties  are  baaed  on  the 
reference  temperature 


(«f  yO*  •  0.664 


Tt*?3 


(Be  -  17) 


U 


1 

2 

3 

4 

5 

6 
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12 
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I* 
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% 
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% 

n 
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99 
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0.690 
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-0.7 
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4J 

101 

6.4 2 

O.651 

0.650 

-0.2 

800 
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-1.0 

0.839 

0.834 

0.6 

2d 

99 

8.17 

0.620 

0.619 

-0.2 

1182 

1155 

-2.5 

0.828 

0.820 

1. 

3a 
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9.57 

0.593 

0.592 

-0.2 

1559 

1506 

-3.0 
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0.810 

-0.4 

Id 
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11.67 

0.555 

0.555 

-0.0 

2246 
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-5.3 

0.807 

O.8O3 

0.5 

6a 
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15.45 

0.501 

0.505 

0.8 

3565 

3322 

-7.0 

0.798 

0.797 

0.1 

7a 

90 

20.45 

0.462 

0.470 

1.7 

5135 

4870 

-5.1 

0.781 

0.793 

-1.5 

3« 

S|00 

0.22 

0.674 

0.664 

-1.5 

400 

403 

41 

402 

2.43 

0.629 

0.630 

0.2 

800 

791 

-2.2 

0.840 

0.827 

1.6 

2a 

400 

3.47 

0.601 

0.600 

-0.2 

1182 

1159 

-3.0 

0.828 

0.815 

1.6 

5b 

401 

4.28 

0.575 

0.576 

0.2 

1556 

1516 

-3.5 

0.813 

0.809 

0.6 

le 

405 

5.52 

0.541 

0.54? 

0.2 

2246 

2154 

-5.1 

0.806 

0.603 

0.4 

6d 

401 

7.49 

0.499 

0.495 

-0.8 

3565 

3366 

-6.3 

0.796 

0.797 

-0.1 

7b 
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9.5  V 

0.456 

0.456 

-0.4 

5135 

4898 

-4.8 

0.780 

0.792 

-.5 

8o 

398 

12.44 

0.419 

0.414 

-0.7 

8160 

7868 

-3.7 

0.786 

9» 

396 

15.95 

0.383 

0.375 

-2.1 

12645 

4e 

794 

0.26 

0.671 

0.664 

-1.0 

800 

2f 

800 

1.74 

0.638 

0.635 

-0.5 

1182 

1190 

2.1 

O.8I3 

0.812 

0.1 

5o 

808 

2.45 

0.612 

0.610 

-0.3 

1550 

1558 

1.1 

0.816 

0.806 

1.2 

if 

797 

3.52 

0.576 

0.574 

-0.3 

2246 

2235 

-1.0 

0.788 

0.803 

-1.9 

6f 

810 

4.95 

0.533 

0.527 

-1.1 

356 5 

3510 

-1.6 

0.806 

0.797 

1.1 

7f 

792 

6.45 

0.490 

0.490 

0.0 

5135 

5162 

0.8 

0.796 

0.792 

0.0 

8a 

805 

8.54 

0.448 

0.443 

-1.1 

8160 

8135 

-0.3 

0.786 

9 1 

76 5 

11.33 

0.405 

0.399 

-1.5 

12645 

Jau«  runs  were  not  fully  evaluated  since  the 

property  values  ware  not  known 

for  th«m. 


Table  la,  Wall  at  reoovery  temperature 
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Run 

’Tr. 

Hd 

15d 

20f 

lOf 

Ilf 

16f 

13f 

18f 

lOe 

11« 

16a» 

17* 

22a 

13e 

18a 

19« 

23o 

lOd 

lid 

16d 

17h 

22d 

13d 

18d 

19* 

23« 


known 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

u 

z 

t 

§ 

h 

T. 

OR 

T* 

oft 

Ma 

m 

1 

0 

(CeiRS)a 

Cr  * 

( — L) 
'2St' 

(Pr-)** 

i 

t 

cf 

toq) 

% 

97 

400 

8.63 

0.667 

0.661 

-0.9 

0.911 

0.795 

-u 

0.934 

-17 

101 

400 

11.84 

0.629 

0.634 

-0.8 

0.862 

0.774 

-11 

0.924 

-19 

94 

400 

19.19 

0.584 

0.583 

-0.2 

0.804 

0.754 

-  6 

0.878 

-16 

101 

1182 

8.93 

0.613 

0.614 

-0.2 

0.873 

0.766 

-14 

0.916 

-19 

102 

1182 

12.37 

0.594 

0.589 

-0.8 

0.814 

0.755 

-  7 

0.876 

-16 

105 

1182 

18.74 

0.553 

0.552 

-0.2 

0.776 

0.746 

-  4 

0.835 

-12 

96 

2246 

13.39 

0.553 

0.547 

-1.1 

0.853 

0.742 

-14 

0.898 

-20 

80 

2246 

22.40 

0.524 

0.526 

0. 

0.799 

0.742 

-  8 

0.842 

-14 

402 

1182 

4.55 

0.591 

0.590 

-0.2 

0.722 

0.758 

5 

0.758 

0 

405 

1182 

6.31 

0.572 

0.572 

0.0 

0.704 

0.753 

8 

0.751 

1 

401 

1182 

9.70 

0.538 

0.534 

-0.7 

0.691 

0.750 

9 

0.734 

2 

404 

1182 

14.83 

0.481 

0.478 

—0.6 

0.686 

0.733 

7 

410 

1182 

18.44 

0.450 

0.446 

-0.9 

0.692 

0.731 

6 

398 

2246 

6.68 

0.531 

0.532 

0.2 

0.709 

0.750 

6 

0.744 

1 

394 

2246 

10.20 

0.503 

0.504 

0.2 

0.688 

0.740 

7 

0.720 

2.6 

394 

2246 

15.50 

0.461 

0.458 

-0.6 

0.692 

416 

2246 

18.70 

0.432 

0.437 

-1.1 

0.695 

0.730 

6 

798 

1182 

3.29 

0.621 

0.617 

-0.6 

0.726 

0.758 

4 

0.755 

0 

805 

1182 

4.55 

0.602 

0.604 

0.3 

0.718 

0.750 

4 

0.756 

—1 

792 

1182 

6.99 

0.568 

0.558 

-1.8 

0.708 

0.750 

6 

0.745 

1 

812 

1182 

10.57 

0.513 

0.510 

—0.6 

0.702 

0.735 

5 

784 

1182 

13.43 

0.476 

0.477 

0.2 

0.677 

0.731 

8 

801 

2246 

4.77 

0.574 

0.559 

-2.6 

0.725 

0.750 

3 

0.755 

-1 

804 

2246 

7.22 

0.536 

0.524 

-2.2 

0.672 

800 

2246 

10.98 

0.488 

0.488 

0.0 

0.694 

0.732 

5 

792 

2246 

13.63 

0.462 

0.464 

0.4 

0.685 

0.730 

6 
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for  thna. 


Tablft  lb.  Coolftd  wall. 
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Die  friction  parameter  baaed  on  stream  conditions  is  then 


(Ba  .  18) 

tfith  the  known  values  of  the  stream  temperature  and  the  reference  temperature 
Equation  Be-18  can  be  used  to  determine  the  friction  parameter  as  based  on  stream 
conditions.  This  parameter  is  compared  in  Tables  la  and  b  with  the  result  of  the 
calculations  by  Young  and  Janssen,  and  the  error  in  per  cent  is  entered  into  the 
sixth  or  seventh  column.  It  is  seen  that  the  maximum  difference  between  the  friction 
factors  calculated  by  both  procedures  exceeds  only  in  three  cases  slightly  the  value 
2%.  Generally  the  error  is  probably  within  the  accuracy  of  the  boundary  layejj" 
solutions,  figure  4  contains  as  dashed  lines  also  friction  parameters  (Cf  ^fR?)s 
which  have  been  calculated  with  incompressible  relationship  Ba  -  1  and  the  reference 
temperature  as  given  by  Equation  Ba  -  16  with  the  same  procedure  as  used  in 
connection  with  Table  1.  It  can  be  recognized  that  the  agreement  between  the  exact 
calculation  and  the  simple  calculation  procedure  using  the  reference  temperature  is 
excellent  regardless  of  the  specific  law  for  the  viscosity  variation. 

The  papers  which  oonsider  the  temperature  variation  of  all  property  values 
report  on  recovery  factors  which  can  not  be  expressed  by  the  simple  Equation  Ba  -  5. 
References  41  and  59  use  the  following  definition  far  the  temperature  recovery  f acton 


r 


0*  -  19) 


figure  6  shows  these  recovery  factors  and  indicates  a  considerable  variation  with 
bach  number  and  stream  temperature.  This  large  variation  of  the  recovery  faster  is 
surprising  and  it  may  be  supposed  to  be  caused  by  the  foot  that  according  to  the 
definition  Equation  Ba  -  19  the  difference  between  recovery  temperature  T  and 
static  temperature  T.  is  compared  with  the  difference  between  total  t«npe$ature  and 
stream  temperature  which  a  fluid  would  have  if  it  had  a  constant  specific  heat  equal 
to  its  real  value  at  stream  condition.  Actually,  the  specific  heat  varies  through¬ 
out  the  boundary  layer  and  an  average  value  should  be  used  to  determine  the 
difference  between  total  and  static  temperature.  It  may  be  calculated  how  well  this 
average  value  for  the  specif lo  heat  is  approximated  by  the  specific  heat  which 
corresponds  to  the  reference  temperature  given  in  Equation  Ba  -  16.  Therefore  the 
recovery  values  in  Figure  3  were  recalculated  (Ref.  23)  using  the  following 
definition  for  the  recovery  faotor 


(Ba  -  20) 
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% 

Mo 

(«f  |G>( 

1182 

4.55 

0.591 

6.31 

o.572 

9.70 

0.538 

14.43 

o.481 

2246 

6.68 

o.531 

10.20 

O.503 

15.50 

0.461 

o.598 

o.3U 

-0.881 

o.738 

0.319 

-0.525 

0.684 

o.325 

-o.4l7 

e.752 

0.313 

-0.321 

1.460 

o.312 

-1.365 

1.085 

o.319 

-O.640 

0.967 

0.305 

-O.450 

These  recovery  factor*  arc  also  plotted  in  Figure  6  end  it  oan  be  seen  that  meat  of 
the  variation  of  the  recovery  factor  with  Mach  number  and  stream  temperature  has 
disappeared.  The  remaining  variation  is  wall  expressed  by  Iquation  Be  -  5  when  the 
Prandtl  number  is  introduced  at  the  reference  temperature  described  by  Squat ion 
Be  -  16.  This  can  be  recognised  from  table  la  in  whioh  the  recovery  temperature  as 
calculated  by  Young  and  Janssen  is  eontainsd  in  column  7  and  the  recovery  temperature 
detennined  with  equations  Ba  -  5  and  Be  -  16  in  column  8.  The  difference  in  per  oent 
is  entered  in  column  9*  It  exceeds  only  in  a  few  cases  5  per  cent. 

Crocoo  already  pointed  out  (Ref.  3k)  that  the  relationships  obtained  for  a  gas 
with  constant  specific  heat  remain  materially  unchanged  for  a  gas  with  variable 
specific  heat  when  the  product  c_dT  is  replaced  by  the  differential  of  the  enthalpy 
i.  In  newer  papers  (Ref.  142  and  36)  this  fact  is  utilised. 

Van  Driest  bases  his  recovery  factors  on  enthalpy  instead  of  temperature  and 
uses  the  following  definitions  for  the  enthalpy  recovery  factor 


r  -  ir  —  ls 
it  -i. 


(Ba  -  21) 


Zt  oan  be  observed  that  this  definition  is  very  similar  to  the  ode  given  in  the 
above  Iquation  Ba  -  20.  Figure  7  contains  as  solid  lines  the  enthalpy  recovery 
factors  as  calculated  in  Reference  36  for  two  different  conditions  in  the  stream 
outside  the  boundary  layer.  The  upper  graph  of  the  figure  presents  the  recovery 
factor  for  a  stream  with  a  statio  temperature  of  400  °R.  The  lower  diagram  gives 
the  recovery  factor  for  a  total  temperature  in  the  stream  equal  100°/.  a  statio 
temperature  of  400°R.  was  chosen  to  approximate  conditions  as  encountered  by  an 
aircraft  in  flight  through  the  atmosphere  whereas  a  total  temperature  of  100°F.  is 
supposed  to  represent  conditions  as  found  in  wind  tunnels.  A  comparison  of  the 
recovery  factors  in  this  figure  taken  from  Ref.  59  with  the  enthalpy  recovery 
factors  in  Table  la  discloses  a  marked  difference  especially  at  higher  Maoh  numbers. 
This  is  due  to  the  different  temperature  variation  of  FTandtl  number  on  which  the 
calculation  by  Young  and  Janssen  on  on*  side  and  by  Van  Driest  on  the  other  side  is 
baaed.  Zt  was  pointed  out  before,  that  the  lower  FTandtl  masher  values  were  used 
in  Ref.  59,  whereas  Van  Driest  adopted  the  set  of  higher  values  in  Figuze  3.  Also 
inserted  in  the  figure  as  dashed  lines  are  recovery  factors  which  have  been 
calculated  uaing  Iquation  Ba  -  5  into  which  the  FTandtl  mniisr  was  introduced  at  a 
reference  temperature  as  given  by  Iquation  Ba  -  16.  Zt  mgr  be  observed  that  the 
agreement  between  the  exact  calculation  and  the  calculation  from  Equation  Ba  -  5  is 
very  satisfactory.  The  maximum  deviation  between  both  sets  of  calculations  is  less 
than  2  per  oent.  ITom  the  data  publiahed  in  Reference  Ba  -  28  the  enthalpy  recovery 
factors  have  been  calculated  inserted  in  Sable  la  as  column  10.  They  oan  be 
compared  with  the  values  9  Fr*  in  column  11.  The  difference  in  per  oent  is  entered 
into  column  12.  Zt  may  be  observed  that  the  difference  is  only  in  one  ease  larger 
than  2  per  oent.  Zt  is  smaller  than  the  error  in  column  6  for  the  temperature 
recovery  factors.  Use  of  the  enthalngreoovery  factor  is  therefore  preferable.  Zt 
is  also  remarkable  that  the  value  approximates  the  actual  recovery  factor  so 

well  regardless  of  the  FTandtl  number  variation  on  which  calculations  have  been 
based. 
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Heat  tranafar  coefficients  have  bean  calculated  by  Young  and  Janssen  for  the 
aaae  range  of  parametera  aa  the  friction  faotora.  One  result  of  this  as  of  all 
other  calculation  procedures  is  that  the  Nueeelt  number  as  a  diaansionlsas  paras* ter 
for  the  heat  transfer  coefficient  increases  proportional  to  the  square  root  of  the 
Reynolds  number.  Acoor^naly,  Young  and  Janssen  (Ref.  59)  presented  their  data  aa 
the  parameter  Nu/  *f!H?  f!P  in  which  all  the  property  values  have  been  Introduced 
at. the  stream  temperatures  assuming  that  the  Nuaselt  number  is  also  proportional  to 


Since  it  has  been  shown  that  the  friction  factor  can  be  expressed  very 
accurately  by  a  simple  formula,  it  appears  advantageous  to  calculate  the  heat 
transfer  coefficient  with  the  help  of  the  friction  factor,  therefore  the  ratio 
«f/2  St  which  according  to  Iquation  Be  -  8  has  the  value  Fr  for  a  fluid  with 
constant  property  values,  waa  determined  frot  Young  and  Janssen's  data,  tfith 
property  values  based  on  stream  temperature  toe  following  aquation  holds 


•  •  a  (ff/Vfe  r.  <»** 

(Be  -  22) 

The  parameters  on  the  right  side  of  the  equation  are  tabulated  in  Reference  59* 

Now  it  will  be  checked  how  well  Iquation  Ba  -  6  approximates  the  values  calculated 
by  Young  and  Janssen  when  all  properties  in  this  equation  are  introduced  at  the 
reference  temperature  described  by  Iquation  Ba  -  16.  With  Iquatlons  B  -  3  and 
B  -  6  the  ratio  friction  factor  to  Stanton  number  can  bo  written 


(Sf-  \  . 

lt;*  h  Te 

(Ba  -  23) 

The  only  property  value  in  this  equation  is  the  specific  heat  0  .  Therefore  the 
corresponding  ratio  based  on  properties  at  reference  temperature  is 


•  (fn>. 

(Ba  -  24) 

This  parameter  aa  calculated  from  Young  and  Janssen's  data  is  comtainsd  in  ooluan  8 
of  Table  lb.  The  values  Fr«t/3  are  contained  in  column  9  and  the  percentage 
difference  in  column  10.  Zt  can  be  observed  that  the  agreement  for  the  heat 
transfer  parameter  is  not  as  good  as  for  the  friction  factors  and  the  temperature 
recovery  factors.  Only  for  the  two  stream  temperatures  4 00°R.  and  800°R.  is  it 
within  10  per  cent.  Calculations  by  Hunker  and  KeoLean  (Ref .  41  and  42)  are  less 
extensive  than  the  one  in  Reference  59  but  are  generally  in  good  agreement  with  the 
latter.  Therefore  a  comparison  with  the  incompressible  relationship  was  not  made 
for  these  calculations. 
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Tan  Driest*  has  also  calculated  beat  transfer  coefficients  with  the  Jfraadtl 
number  ▼ariatioa  as  proposed  by  the  National  Bureau  of  Standards.  Ha  used  in  the 
presentation  of  his  results  a  Stanton  lumber  which  is  based  on  enthalpies  i  instead 
of  t superstores  analogously  as  in  the  definition  of  the  recovery  factor.  Instead 
of  the  customary  definition  of  the  Stanton  number 


St 


t  ®P  *e  (*r  “  *s) 


he  defines  it  in  the  following  way 


(Ba  -  25) 


at.  .  — -■  % 

(Ba  -  25m) 

She  paper  contains  plots  of  the  parameter  c«/2  St  whi^afer  oonstant  Prandtl 
number  has  according  to  fetation  Ba  -  8  the  value1  (ft: )“  .  figure  8  shows  this 
parameter  plotted  over  the  Mach  number  for  two  different  conditions  in  the  free 
stream.  As  upper  diagram  holda  for  a  free  stream  static  tnaparnlurs  of  BOO*.  and 
the  lower  diagram  holds  for  a  total  free  stream  temperature  of  100°JT.  *  560*. 
Compared  with  these  boundary  layer  solutions  are  again  the  resulte  of  am  spgromftmate 
procedure  shioh  sets  the  psraamtar  e«/8  St.  equal  to  Hr*  and  Introduces  the  Hraadtl 
aaahsr  at  the  refsrease  temperature  am  glean  w  Bjuetlom  Bt  -  16.  It  any  he 
observed  that  the  — rooaont  bstwsam  both  sets  of  calculations  is  again  amcelleat. 

The  deviation  is  in  the  whole  range  investigated  less  thorn  1  par  sent.  In  view  of 
this  good  agreement  tho  par  asm  ter  f/2  St.  was  also  ealeulated  from  Young  and 
Jbaooon's  data.  It  is  oontaiasd  in  aolwsi  7  of  fable  lb.  flnlunn  8  gives  again  tbs 
par  cent  difference  between  the  actual  values  in  coluna  7  and  the  appreagmatioa  by 
Ipiatioa  Ba  .  8  as  contained  in  mil— nn  $.  ft  earn  bo  obenrrad  that  the  opaaaoat  is 
poor  for  a  stream  temperature  T.  -  100*.  and  800*.  A  stream  tanmintiue  of  100*. 
nay  occur  la  wind  tunnel  toots  mat  is  asvey  ana  numbered  in  the  atanaphase.  Be 
design  calouletioas,  therefore,  the  use  of  the  Stanton  nanbrnr  based  on  enthalpy  is 
preferable.  Bis  Stanton  matter  can  be  determined  from  fetation  Ba  -  8  with  on 
error  less  than  2.8  par  seat,  then  calculations  of  rsooesry  fastoro  or  heat 
treaafsr  ooafflolaata  are  baaed  oa  enthalpy,  it  is  probably  mors  loaoealsat  to 

a  refsrease  enthalpy  instead  of  tho  refanaass  tsnpsraturs  (feist  I  on  *  -  14).  Bar 
this  roanon  the  oaloalation  prtiaedais  as  is  scribed  in  gpiallr  Z  man  repeated  on 
the  basis  of  anbhalpdas  instead  of  tsapsreturss.  It  was  found  that  within  thn 
accuracy  obtainable  the  equation 


i*  -  i,  ♦  o.5  (i*  -  1J  ♦  e.2t  (1  -  i» 

^  ^  <*C  -  *) 


*) 


p.ia5. 


30 


which  hu  the  im  numerical  constant*  as  Iquation  Ba  -  16,  datermines  tha  rafaraaoa 
conditions  whioh,  whan  os  ad  for  tha  propartiaa  in  tha  aquation  for  tha  friction 
factor,  gives  bast  agreement  with  exact  boundary  layer  solutions.  The  refer* noa 
enthalpy  nay  also  be  expressed  in  a  different  way  ass 

i#  -  ».5  (i,  ♦  V  ♦  0.22  r±t£  (Ha)2  i, 

(Ba  -  26a) 

A.  comparison  of  measured  temperature  recovery  factors  and  heat  transfer 
ooeffloients  with  calculated  values  has  been  made  in  a  recent  paper  by  0.  Ibar 
(Bef.  155).  Sber  shows  that  measured  temperature  recovery  factors  in  a  laminar 
boundary  layer  on  cones  vary  around  the  value  0.Q5  with  a  maximum  deviation  of 
1  per  cent.  Within  this  limit  the  values  agree  also  with  the  calculated  results 
which  have  been  presented  in  this  paper. 

In  the  Figure  k  of  his  paper.  Bber  compared  measure  recovery  factors  with 
values  whioh  have  been  calculated  by  Hunker  and  Ms  La an  and  points  out  that  for 
larger  Mach  numbers  a  considerable  difference  exists.  Zt  has  to  be  kept  in  mind, 
however,  that  the  calculated  data  by  Hunker  and  McLean  have  beau  obtained  for 
free-streaas  conditions  which  were  different  frost  the  ones  for  whioh  the  measure¬ 
ments  have  been  seeds,  The  measured  recovery  factors  have  to  be  compered  with  the 
values  in  Figure  7  and  exhibit  then  an  mpeamcat  which  again  is  within  1  per  oast. 
The  peculiar  fact  that  temperature  recovery  factors  msaaured  on  a  flat  plate  were 
consistently  higher  than  the  ones  determined  tram  comes  and  aenmnrt  values  around 
0.86  has  been  explained  by  Korobkin  (Bef.  Ji3)  as  caused  by  heat  conduction  in  the 
solid  plate  material.  Heat  transfer  coefficients  could  only  he  ammmirod  with  an 
accuracy  of  approximately  lo  per  cent  and  within  this  limit  again  ^pmamsnt 
between  values  measured  on  cones  and  calculated  heat  transfer  coefficients  have 
been  observed. 

Beconmeaded  Calculation  lbocedure 

On  the  basis  of  the  investigation  made  in  this  chapter  the  following  isnneibns 
ia  recommended  for  the  calculation  of  the  friction  and  heat  transfer  of  flat  plates 
in  supersonic  laminar  flows 

The  wall  shearing  stress  ia  calculated  from  the  following  equation 


(Ba  -  27) 

with  tha  friction  factor  appearing  in  this  equation  given  by  the  relationship 


cf 


(Ba  -  26) 
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The  property  values  are  introduced  into  both  of  the  above  equations  at  a  reference 
temperature  which  may  be  tentatively  calculated  from  the  following  relationship 


T*  =  Tg+  o.5  (Tu  -  Ts)  «•  o.22  (Tr  -  Tg) 


(Ba  -  2?) 


This  procedure  was  shown  to  give  agreement  with  exact  boundary  layer  aolutions 
within  £  2  per  cent. 

Heat  transfer  coefficients  are  defined  by  the  following  equations  i 


qu  .  h  (Tr  -  Tw) 

for  constant  specif lo  heat,  ori 

<»«  •  hi  (1r  -  ‘v> 

for  variable  specific  heat. 

The  recovery  temperature  is  determined  by  the  relation 


(Ba  -  30) 


(Ba  -  3D 


T  ■  T  ♦ 


,1 

2  cp 


for  constant  speoific  heat  and  by 


(Ba  -  32) 


i 


r 


i.  ♦  r4 


(Bh  -  33) 


far  variable  specific  heat. 

The  recovery  factor  appearing  in  both  of  the  above  equations  is  calculated 
from  the  equation 


r  = 


(Ba  -  3k) 


into  which  the  Jtaadtl  number  has  again  to  be  introduced  at  the  reference  temperature 
as  given  by  Squat  ion  Ba  -  24.  Alternatively  the  irandtl  number  to  be  introduced  into 
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Equation  Ba  -  34  oaa  also  ba  determined  Aron  the  rafaranoa  enthalpy 


i#  •  o.5  (is4»  iw)  0.22  rA  Tzk  (Ma)  i8 

(Ba  -  as) 

Recovery  factors  oalculat* i  in  thia  way  agree  within  ♦  2  par  cant  with  the  results 
of  exact  boundary  layer  solutions. 

1!he  heat  transfer  coefficient  appearing  in  Equation  Ba  -  JO  and  Ba  -  31 
calculated  from  the  Stanton  number  value 


St  *  - h _ 

9  Cp  Va 


for  constant  speoifio  heat  or  fron 


(Ba  -  36) 


(Ba  -  37) 


for  variable  specific  heat.  The  Stanton  nuabers  are  obtainable  fron  the  relationship 


St  »  £  (Pr-j* 

(*  -  3») 

into  which  the  irandtl  numbers  have  again  to  be  introduced  at  the  reference 
tcnperature  given  by  Squat  ion  Ba  -  28  or  the  reference  enthalpy  Equation  Ba  -  35. 
Heat  tranafer  coefficients  calculated  in  this  way  agree  within  2.o  per  cent  with 
the  results  of  exact  boundary  layer  solutions,  for  free  atrean  tanperaturea  above 
300°*. 

The  nain  uncertainty  in  the  calculation  of  recovery  factors  and  heat  trenafer 
coefficient  a  ia  caused  by  our  insufficient  present  day  knowledge  of  Araadtl 
limbers  for  air  at  hijmr  tanperaturea.  irandtl  nunber  values  aa  given  in  the 
Qea  Tables  and  reoo— ended  by  the  National  Bureau  of  Standards  differ  for  instance 
at  a  temperature  of  2000°R.  by  12  per  cent  and  accordingly  the  recovery  f as  tore 
calculated  with  those  Brandt!  numbers  differ  by  6  per  cent  and  heat  transfer 
coefficients  by  8  per  cent. 
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TURBULBVT  BOUNDARY  LAUR  PLOW 


In  tbs  proceeding  chapter  it  has  bean  damonstratad  that  tha  oaloulation  of 
friction  and  ha  at  tranafar  in  laminar  boundary  layar  flow  haa  baan  rafinad  to  an 
ao curacy  which  ia  vary  difficult  to  parallel  in  experimental  investigations.  In 
turbulent  boundary  layers  tha  aituation  la  completely  different.  Tha  present  day 
understanding  of  turbulent  exchange  of  moawntum  and  energy  especially  under 
conditions  of  widely  varying  temperature  ia  restricted  to  such  an  extent  that  exact 
solutions  of  the  boundary  layer  equations  can  not  be  obtained.  In  any  oaloulation 
procedure  assumptions  have  to  be  introduced  on  the  exchange  neohaniam  and  these 
assumptions  have  to  be  continuously  checked  with  experiments.  In  obtaining  friction 
and  heat  transfer  data  we  have,  therefore,  to  rely  heavily  on  measurements. 
Calculations  are  neoessary  and  useful  in  order  to  extend  the  renge  of  pareamters 
outside  the  experimentally  investigated  one.  However,  it  will  be  neoessary  to 
select  from  the  different  calculations  procedures  published  in  the  literature  those 
which  agree  well  with  experimental  data  in  the  investigated  range.  Gonparatlvsly 
extensive  experimental  investigations  are  available  for  the  friotion  factor  of 
plates  whose  temperatures  are  equal  to  the  reoovexy  temperatures. 

In  Figure  9  (Ref.  65)  the  information  on  the  influence  of  Maoh  number  on  this 
friction  factor  is  collected.  The  measurements,  tha  results  of  which  are  presented 
in  the  figure,  have  been  made  on  fist  plates  or  cylinders.  Plotted  is  the  ratio  of 
the  friction  factor  in  compress! hie  flow  to  the  friction  factor  in  incompressible 
(means  constant  property  value)  flow.  The  friction  factor  is  again  based  on  the 
property  values  taken  at  stream  tmeperature .  It  may  be  observed  that  it  decreases 
considerably  with  increasing  Mach  number.  The  influence  of  Reynolds  number  on  the 
friotion  factor  has  been  found  to  be  the  same  for  supersonic  as  for  low  velocity 
flow.  Therefore  Figure  9  applies  to  aay  value  of  the  Reynolds  number,  friotion 
factors  which  have  been  measured  on  cones  by  Rradfleld  (Ref.  18)  have  values  which 
arm  approximately  20  per  cent  higher  than  the  results  of  the  measurements  collected 
in  Figure  9.  The  reason  for  this  discrepancy  is  not  explained  until  now. 

Figure  10  (Ref.  65)  presents  the  predictions  which  different  published  theories 
make  on  the  change  of  friotion  factors  with  Meoh  number.  It  may  be  observed  that 
the  various  predictions  differ  to  a  vary  high  degree.  The  best  agreement  to  the 
experimental  results  presented  in  Figure  9  give  calculations  made  by  Wilson, 
Monaghan,  Oops,  fraakl-Yolahel,  and  Thsher.  lit  even  with  these  results  it  is 
difficult  to  make  predictions  on  the  friction  factor  at  Mmh  numhars  because 

the  different  theories  deviate  considerably  at  Mach  numbers  shame  five. 

Raced  with  this  situation  it  is  appropriate  to  look  for  e  simple  calculation 
proosdure  which  agrees  with  the  experimental  results  in  tha  investigated  range. 

Toung  and  Jeneeea  (Rtf.  59)  pointed  out  that  the  noo  of  constant  property  relations 
with  tha  refsrenss  temperature  established  for  laminar  boundary  layers  gives  good 
agreement  for  turbulent  boundary  layers  as  wall.  War  this  procedure  the  relation¬ 
ship  for  akin  friotion  in  n  constant  property  fluid  flow  must  be  eatahliahed  at 
first.  Several  equations  are  offered  in  the  literature  for  this  purpose. 

The  HLaalus  aquation 


~  *  0.0296  Re"0,2 
2 
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(*>  -  1) 


gives  good  agreement  with  measured  values  up  to  Ra  »  10?.  At  larger  Reynolds 
□umber?  Situation  Bb  -  1  gives  friction  factors  which  are  too  low. 

A  good  representation  of  tha  measured  average  friction  factors  in  the  whole 
range  up  to  Re  *  10?  ia  obtained  by  the  Karman-Sehoenherr  equation 


-  ioglQ  (r«  ) 

Y  Cf  (Bb  -  2) 


The  local  friction  faotor  may  be  calculated  from  the  average  value  by  the 
relation 


°f 


q.557 
0.557  2 


(Bb  -  3) 

Simpler  to  apply  is  a  formula  by  Sohuls-CTunow  for  the  local  friction  faotor 


cf  - 


o«370 


(logi0  He) 


2.584 


(Bb  -  4) 

and  a  relation  by  Praodtl-Sohlichting  for  the  average  friction  faotor 


cf 


0.455 

- - “735 

U°gl0  Re) 


Both  of  these  relations  represent  measurements  up  to  Be 


(Bb  -  5) 
IQ?  vary  wall. 


Since  the  measurements  plotted  in  Figure  9  have  been  made  at  Beynolds  numbers 
below  10' .  the  BLaslus  Equation  (Bb  -  l)  will  be  used  to  establish  the  Nash  number 
dependence  of  the  friction  faotor  with  the  use  of  the  reference  temperature-  The 
friction  footers  in  Figure  9  are  baaed  on  property  values  at  stream  temperature 


(Bb  -  6) 
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and  compared  with  tha  friction  factor  ofi  of  an  incoopresaible  (constant  property) 
fluid  at  stream  temperature 


cf  =  °-0296 

2  1  <Ss  V9  *  /^8)0,2 

(8b  -  7) 

Assuming  that  for  a  variable  property  (compressible)  fluid  the  sheering  stress  is 
given  by  the  constant  property  value  relationship  when  the  property  values  are 
Introduced  at  the  reference  temperature  T*  results  in  the  equation 


rv  = 


o . 0296 

Vs  x  //> 


oT2 


(Bb  -  8) 


A  comparison  between  the  tquations  8b  -  6,  fib  -  7.  end  Bb  -  8  gives  the  following 
formula  by  which  the  ratio  of  the  friction  faotor  in  supersonic  flow  to  the  friction 
factor  for  a  fluid  with  constant  property  value  can  be  calculated. 


(8b  -  9) 

The  results  of  such  a  calculation  in  which  it  was  aasu—ri  that  the  viscosity  varies 
promotional  to  the  power  of  0.76  of  the  absolute  teaperature  is  indicated  by  • 
dashed  line  in  figure  9.  Zt  nay  be  observed  that  this  dashed  line  average*  the 
measured  values  very  well.  It  also  agrees  at  higw  Mach  numbers  with  the  calcula¬ 
tions  fey  Wilson,  ftonaghsa,  and  ftioker  as  presented  in  figure  10.  Zt  can  therefore 
be  recaanwnded  for  a  plausible  extrapolation  of  the  test  results  into  the  higher 
Mach  number  range.  A  pap**r  (by  I.  H.  Abbot  -  Ref.  181)  presented  to  the  fourth 
meeting  of  A0AB8  wind  tunnel  and  model  testing  panel  oontains  the  results  of  test 
on  thin-walled  tubes  in  a  firing  rang*.  TWo  series  of  test  points  ware  obtained, 
one  at  Maoh  number  3.9  and  for  the  wall  cooled  to  stream  static  temperature,  the 
other  at  Itooh  number  7*25  and  3L,  ■  1.8  I  .  loth  groups  of  tost  points  are  abown  in 
figure  9.  The  deshad  lias  in  this  figure  represents  tha  ratio  *■  ealaulatad 

with  ^rations  Mb  -  9  and  la  -  16.  The  agreement  la  especially  goodror  tha  meal  Tar 
Nash  number,  for  which  according  to  Abbott  tha  tests  reproduced  moot  satisfactorily. 

Scan  friction  coefficients  on  a  cooled  well  are  also  report ud  in  Refers noe  180. 
They  were  measured  at  Hash  numbers  between  3  end  7.7.  The  values  Of/ofi  for  these 
tests  are  all  lower  than  even  the  curve  for  the  wall  at  recovery  t  separators  shown 
in  Figure  9.  Zt  la  believed  that  this  is  caused  tgr  tha  fast  that  tha  measurements 
ware  made  on  a  wind  tunnel  wall  along  which  tha  pressure  was  not  constant  but 
decreased  in  flow  direction. 
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RifwtDM  73  contains  the  reeulte  of  e  very  extensive  investigation  which  has 
been  ande  at  the  NiCA  to  study  friction  factors  and  heat  transfer  ooefflclente  for 
air  flowing  turbulestly  through  a  round  tube  with  large  teaperature  difference*  and 
eubeonio  velocities.  The  essential  result  of  these  tests  was  that  the  normal  oonatant 
property  value  relationship  described  the  results  of  these  tests  satisfactorily  whan 
the  property  values  were  introduced  at  a  reference  teagerature  which  was  situated 
half  wqr  between  the  wall  teaperature  aad  the  air  teaparature.  Zt  aay  be  observed 
that  equation  Be  -  16  contains  the  seme  relation  for  the  reference  teaperature  for 
the  case  that  the  sir  velocity  is  small  (recovery  teaperature  equal  to  wall 
temperature).  Therefore,  it  aqr  be  Justified  to  reooammd  use  of  oonstant  property 
value  relationships  togsther  with  equation  ha  -  16  as  reference  temperature  for 
turbulent  flow  along  a  flat  plate. 

Ifer  plausible  deductions  Aokemann  (Ref.  62)  and  Squire  (Ref.  83)  determined 
the  following  relationship  for  the  recovery  factor  in  a  turbulent  boundary  layer 


r  -  V*r 

(Bb  -  10) 

Test  reeulte  on  flat  plates  as  well  as  on  cones  agree  very  well  with  this 
prediction  so  that  it  is  today  in  eoamon  usage  to  describe  the  teaparature  recovery 
in  a  turbulent  boundary  layer  on  a  flat  plate.  The  test  results  showed  no  conclusive 
variations  of  the  recovery  factor  with  Reynolds  number  or  with  Mach  number.  A 
slight  decrease  of  the  reoovery  factor  with  Mach  number  which  has  been  predicted  ty 
an  extension  of  Squire's  calculation  into  the  compressible  flow  rang*  (Ref.  87)  has 
up  to  uow  not  been  confirmed  by  '’xp«riments,  probably,  because  this  effect  is  too 
snail  in  the  range  0f  Mach  numbers  investigated  till  now.  It  has  to  be  remembered 
that  all  experimental  information  has  been  obtained  at  a  total  temperature  of  the 
free  stream  around  100°F. 

The  simplest  way  to  calculate  heat  transfer  in  turbulent  flow  is  again  to 
derive  it  from  the  friotion  factor  with  an  equation 


analogous  to  relationship  Be  -  8  in  which  the  function  S  has  to  be  determined 
experimentally . 

Jbr  email  tamper ature  differences,  end  low  velocities  Colburn  (Ref.  66)  has 
found  that  this  factor  S  can  be  expressed  for  turbulent  flow  as  well  ee  for 
flow  as  ljtr  *'3  Accurate  test  results  for  which  this  expression  could  be 

checked  in  high  velocity  flow  with  large  temperature  differences  are  not  known  to 
the  author. 
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Rubeein  (Ref .  82)  obtained  in  a  oaloulation  in  which  ha  had  again  to  aaka 
plauaibie  assumptions  on  the  turbulent  exchange  aechanieae  the  raault  that  the 
factor  S  agrees  for  high  velocity  high  temperature  difference  flow  for  a  Reynolds 
number  below  10®  with  the  expression  1 /St**  and  that  it  deoreaeea  elightly  with 
Raynolde  number  ao  that  for  a  Raynolde  number  of  107  it  aeeumes  for  air  the  numerical 
value  1.20  and  for  a  Raynolde  number  equal  109  the  value  1.18.  The  influence  of 
Mach  number  and  of  the  wall  to  atream  temperature  ratio  wee  in  thia  oaloulation 
found  to  be  extremely  email. 

Van  Drieat  (Ref.  36)  makes  in  his  well  known  oaloulation  of  heat  transfer  in 
the  turbulent  high  velooity  boundary  layer  the  assumption  that  the  Arandtl  number 
and  the  ratio  of  turbulent  diffusivlties  for  momentum  and  heat  are  both  equal  to 
one.  for  such  a  fluid  Reynolds  analogy  is  exactly  fulfilled.  This  means  that  for 
Me  ■  4  the  frlotion  factor  of  a  wall  with  the  temperature  equal  to  the  recovery 
temperature  as  calculated  by  Van  meat  is  approximately  20  per  cent  larger  than  the 
one  measured  experimentally.  Zt  is  therefore  understandable  that  the  heat  transfer 
coefficient  calculated  by  Van  meet  agrees  well  with  experiments  and  with  the 
predictions  of  Rube sin. 

ill  the  discussion  in  this  chapter  assumed  that  the  plate  surface  is  smooth. 

Zt  is  known  that  roughness  increases  friction  and  heat  transfer  in  turbulent  flow 
when  the  roughness  exceeds  certain  limits,  extensive  systematic  Information  on 
rough  surfaces  exists  only  for  the  frlotion  coefficient  in  low  velooity  flow.  1 
good  survey  over  this  field  is  for  instance  contained  in  Reference  27. 

Suamsrlsing  our  discussion  on  turbulent  boundary  layers  the  following 
recommendation  la  made  for  a  calculation  procedure  t  Constant  property  relationships 
should  be  used  to  calculate  the  friction  factor,  the  recovery  factor,  and  the  heat 
transfer  coefficient,  namely  Equations  Bb  -  1  to  fib  -  5  for  the  friction  factor, 
Equation  Bb  -  10  for  the  recovery  factor,  and  Equation  Bb  -  11  for  the  heat  transfer 
coefficient.  The  factor  S  in  the  latter  equation  varies  for  air  between  1.23  and 
1.18  in  the  Reynolds  number  range  from  105  to  10°.  Property  values  should  be 
introduced  into  theee  equations  at  the  reference  temperature  given  by  equation 
Ba  -  16  and  the  recovery  factor  as  well  as  the  Stanton  number  describing  the  heat 
transfer  should  be  baaed  on  enthalpies  instead  of  temperatures  when  the  temperature 
range  is  such  that  the  epeoifio  heat  varies  considerably  within  the  boundary  layer. 
The  reooamndatlons  have  to  be  considered  as  tentative  and  have  to  be  oheoked.  by 
more  experiments  at  large  Mach  numbers  and  especially  at  large  temperature 
differences  and  high  temperatures. 
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SECTION  17 


INFLUENCE  07  IRES  SURE  AND  TEMTOA3VRE  VARIATION  ALOMJ  THE 
SURFACE  ON  HEAT  TRAN3FER 


VARIABLR  PRRSSfJRE 


Up  to  now  the  pressure  was  assumed  to  be  constant  along  the  surfaoe  of  the  body 
to  which  heat  is  transferred.  Actually,  the  pressure  usually  varies  along  the  body 
as  well  as  along  the  airfoils  of  aircraft  and  it  is  important  to  know  the  influence 
of  such  a  pressure  variation  on  heat  transfer. 

Exact  solutions  of  the  boundary  layer  equations  for  laminar  flow  and  for 
fluids  with  constant  properties  exist  for  certain  types  of  veloolty  variation 
especially  for  a  velocity  outside  the  boundary  layer  which  increases  proportional 
to  some  power  of  the  distance  from  the  leading  edge. 

Approximate  procedures  which  use  the  integrated  momentum  and  energy  equation 
have  also  been  worked  out  for  the  oaloul  ution  of  heat  transfer  on  surfaces  along 
which  the  pressure  varies  in  some  arbitrary  way.  Information  obtainable  from 
calculations  is  ouch  more  restricted  for  high  velocity  flow  in  fluide  with  variable 
properties.  Only  in  a  very  few  eases  solutions  of  the  boundary  layer  equations  have 
been  obtained  and  even  procedures  which  use  the  integrated  boundary  layer  equations 
become  rather  tedious. 

Generally  it  is  known  that  a  pressure  variation  influences  heat  transfer 
coefficients  to  a  lower  degree  than  the  skin  friction  parametere.  This  is  eaused 
by  the  feet  that  the  pressure  gradient  enters  the  energy  equation  only  in  an 
indirect  way.  In  addition  to  this  fact,  the  shape  a,  used  for  aircraft  bodies  and 
airfoils  designed  to  fly  at  high  velocities,  have  to  be  made  very  stresadined  in 
order  to  reduce  the  drag.  On  such  configurations  which  for  supersonic  flow  have 
also  sharp  noses  or  leading  edges  the  pressure  variation  is  comparatively  small  and 
therefore  the  effect  on  heat  transfer  is  of  minor  importance.  This  finding  is 
substantiated  by  experiments  contained  in  Reference  65  in  which  averts  heat  transfer 
coefficients  were  measured  on  a  cylindrical  body  with  a  conical  nose.  Xhe  pressure 
along  the  cylinder  could  be  varied  by  an  adjustment  of  the  wind  tunnel  walls.  In 
addition  to  a  constant  pressure  along  the  surfaoe.  pressure  variations  were 
investigated  corresponding  to  Maoh  number  variations  along  the  surfaoe  ae  shown  in 
Figure  11.  It  is  mentioned  in  the  referenced  report  that  such  a  variation 
corresponds  approximately  to  the  pressure  variation  as  found  around  a  3  per  cent 
thickness  beconvex  air-foil  or  to  the  pressure  variation  around  a  parabolic  are 
body  of  revolution  with  a  length  to  diameter  ratio  equal  to  15*  It  was  found  that 
for  pressure  variations  as  indicated  in  Figure  11  the  average  heat  transfer 
coefficient  varied  only  to  a  degree  which  wee  analler  than  the  accuracy  with  which 
heat  transfer  coefficients  could  be  measured. 

The  same  conclusion  was  reached  in  Reference  26  where  local  heat  transfer 
values  and  recovery  factors  have  been  aeaeured  on  cones  with  straight  and  parabolic 
generatrices.  It  can  therefore  be  expected  that  for  bodies  with  ill  angle  noses 
or  leading  edges  end  for  flow  under  mil  angle  of  attack  the  ooostant  pressure 
relatione  give  heat  transfer  date  of  sufficient  accuracy  for  desist  purposes.  VI th 
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r aspect  to  cooling  a  rounded  nose  would  be  advantageous  for  missiles  sinoe  in  this 
way  the  large  heat  transfer  coefficients  near  leading  edges  are  avoided.  Heat 
transfer  on  such  a  rounded  nose  will  deoidedly  be  different  from  that  on  a  flat 
plate.  An  extensive  treatment  of  the  influence  of  pressure  variation  on  heat 
transfer  to  bodies  of  arbitrary  shape  would  exceed  the  aoope  of  this  paper,  fbr 
this  reason  no  further  discussion  on  the  influence  of  pressure  variations  will  be 
mode  and  the  reader  is  referred  to  the  literature  presented  in  the  appendix. 

One  remark  may  be  added,  namely  that  for  flow  with  pressure  gradients  the 
correlation  of  the  actual  friction  and  heat  transfer  parameters  with  the  oonstant 
property  relations  ty  speoifioation  of  a  reference  temperature  is  not  posaible.  This 

can  be  se*»n  from  H®fer«nce  99  where  calculations  are  made  under  the  assumption 
that  the  product  of  density  and  viscosity  is  independent  of  temperature.  Friction 
factor  and  Musselt  number  are  found  to  vary  with  the  ratio  wall  to  stream  tempera¬ 
ture.  "o  such  variation  would  be  indicated  by  tb»  constant  property  relation  for 
any  choice  Qf  +he  reference  temperature. 


VARIABLE  WALL  T^MP'RATUR1! 

That  a  tamper at urs  variation  along  the  surface  of  a  heated  object  is  a 
parameter  of  importance  for  the  heat  tranafar  coefficient  has  been  only  recently 
generally  recognised .  fags  and  Falkner  (Ref.  115)  presented  *or  laminar  flow  an 
approximate  treatment  of  this  problem  as  early  as  in  1931  but  it  was  only  after 
1950  that  intensive  work  was  started  to  investigate  the  dependence  of  heat  transfer 
on  this  parameter.  All  calculations  on  hast  transfer  at  variable  wall  temperature 
ore  based  on  the  fact  that  for  constant  property  values  the  energy  aquation  which 
describes  the  heat  flow  through  e  boundary  layer  into  a  wall  is  linear  and  that 
therefore  general  eolutiona  of  this  equation  can  be  obtained  ty  a  superposition  of 
individual  special  solutions.  This  linearity  of  the  energy  aquation  is  in  laminar 
flow  essentially  still  preserved  when  the  property  values  vary  in  auoh  e  way  that 
the  product  of  density  tinea  viscosity  is  oonstant.  Under  this  condition  and  for  e 
specific  type  of  temperature  variation  along  the  surface  1 


Tw  -  TP  -  C 


(Cb  -  1) 


Chapman  and  Rubeain  (Ref.  113)  obtained  exact  solutions  of  the  laminar  boundary 
layer  equations  describing  the  flow  end  heat  transfer  to  a  flat  plats,  figure  12 
contains  heat  transfer  coefficients  based  on  the  local  teonarature  difference 

as  calculated  in  this  reference.  The  parameter  h/VS?  is  plotted  over  the 
exponent  which  determines  the  temperature  variation  along  the  surface.  It  any  be 
observed  that  with  increasing  value  off  a  considerable  increase  in  the  local  heat 
tranafar  coefficient  ns  determined  ty  the  parameter  MU/  VW  occurs.  Mien  the 
temperature  difference  between  the  well  and  recovery  temperature,  which  determines 
the  beet  transfer.  Increases  linearly  with  distance  from  the  leading  edge  the  heat 
transfer  coefficient  will  be  ty  65  per  cent  larger  than  the  heat  tranafar  coefficient 
under  corresponding  conditions  but  for  a  constant  wall  temperature,  fbr  a  temperature 
difference  which  increases  proportional  to  the  square  of  the  distance  from  the 
leading  edge,  the  heat  transfer  coefficient  la  larger  ty  a  factor  of  2  than  the  best 
transfer  coefficient  for  constant  wall  temperature . 
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Conditions  similar  to  the  ones  just  diseusssd  occur  nssr  tbs  boss  or  loading 
sdgs  of  a  body  or  airfoil  of  an  aircraft  in  high  spaed  flight  whan  the  aircraft 
surface  is  cooled  either  ty  radiation  or  by  heat  storage  in  the  solid  material.  In 
this  oase  the  wall  surface  temperature  is  determined  fey  a  balance  between  the  heat 
transferred  from  the  air  to  the  surface  fey  oonrection  and  the  heat  rumored  from  the 
surface  by  radiation  or  fey  heat  storage  in  the  interior.  It  is  known  that  the  heat 
transfer  coefficient  on  the  leading  edge  of  objects  is  rery  large  (theoretically  on 
the  flat  plate  it  is  equal  to  infinity)  and  that  it  decreases  with  increasing 
distanoe  from  the  leading  edge.  Correspondingly,  the  temperature  of  the  object  near 
the  leading  edge  is  practically  equal  to  the  recovery  temperature  and  it  drops  with 
increasing  distanoe  because  the  heat  removal  through  radiation  beoomes  more  and 
more  effective  as  compared  with  the  eonveotive  heat  tranafer.  Temperature  distri¬ 
butions  in  wings  of  diamond  shaped  aa  calculated  fey  bye  (Ref.  116)  show  this 
characteristic  behavior.  It  has  to  be  expected  thi  the  heat  transfer  coefficients 
in  a  close  region  downstream  of  the  leading  edge  of  such  objeota  will  be  considerably 
larger  than  the  heat  transfer  coefficients  calculated  under  the  assumption  of  a 
constant  wall  temperature .  Correspondingly  the  surface  temperatures  will  be  higher. 

It  is  possible  to  determine  heat  tranafer  for  laminar  flow  conditions  and  for 
an  arbitrary  variation  of  the  wall  temperature  by  superposition  of  individual 
solutions  which  have  been  obtained  for  r  temperature  variation  according  to  Iquation 
Cb  -  1  fey  expressing  the  temperature  variation  along  the  surface  as  a  aeries: 


%-*r  ■B’l** 

(Cb  -  2) 

This  method  which  has  been  used  by  different  authors  has  the  disadvantage  that  the 
partial  aolutions  are  known  only  for  laminar  flow  oooditiona  and  for  a  limited 
number  of  terms  in  equation  Cb  -  2.  It  is  often  not  possible  to  empress  a  tempera¬ 
ture  variation  oo curing  on  am  object  fey  a  power  scries  with  this  limited  number  of 
toms  satisfactorily.  Jbr  this  reason  another  method  will  be  described  harm  which 
can  be  applied  to  any  temperature  variation  and  to  turbulent  as  well  as  laminar  flow. 
It  is  mainly  based  on  work  by  Rond  (Ref.  110)  and  Rubesln  (Ref.  120).  this  method 
starts  from  a  solution  of  the  boundary  layer  equations  for  a  step-wise  varying  wall 
temperature  and  replaces  an  arbitrary  wall  temperature  variation  hr  on  infinite 
number  of  infinitely  small  steps,  in  integration  procedure  of  the  known  solution 
far  one  step  then  gives  the  answer  for  the  arbitrary  varying  wall  temperature.  lhla 
derivation  obtains  the  following  equation: 


q(x) 


a^(|,x)^il2 


+  [Tw(  %♦)  -  Tw(*p] 

(Cb  -  3) 


far  *he  heat  flow  at  the  distanoe  X  from  the  lending  edge  into  a  surface.  In  this 
equation  h  is  the  heat  transfer  coefficient  whioh  is  found  at  the  location  X  when 
the  temperature  over  the  surface  is  at  first  equal  to  the  recovery  temperature  and 
then  jumps  suddenly  at  the  location  j  to  a  oertain  value  which  afterwords  is  again 
constant.  T_  is  the  temperature  variation  actually  found  along  the  surface.  The 
integration  Is  to  be  performed  disregarding  discontinuities  in  the  wall  temperature 
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T*.  lb*  summation  expressed  by  tbe  second  term  on  tne  right  hand  side  of  the 
equation  has  to  toe  considered  when  the  wa^.1  temperature  Tw  varies  discontinously 
along  the  surface.  Zn  this  case  (  §j[  )  is  the^wall  temperature  immediately 

downstream  of  a  sudden  temperature  Jump  and  (  §£  )  is  the  corresponding 

temperature  immediately  upstream  of  the  jump.  If  the  temperature  of  the  surface 
differs  at  the  leading  edge  from  the  reoovary  temperature  then  this  summation  has 
also  to  toe  applied  to  the  leading  edge  with  a  temperature  T  (  CT  )  equal  to  the 
recovery  temperature.  Bie  procedure  will  toe  explained  onoe  more  for  a  specif lo 
temperature  variation  later  on. 

The  following  expressions  for  the  heat  transfer  coefficients  which  are  needed 
in  the  above  equation  have  been  found  by  approximate  calculations  for  the  conditions 
on  a  flat  plate,  for  a  laminar  boundary  layer  (Ref.  121) 


h(S,x)  *  VR?  VpP  [l  -  <$)*]"* 

(Cb  -  4) 

for  a  turbulent  boundary  layer  (Ref.  122) 

h(V*>  *  Sigj&Ji  Pr"[l  -  (if]"'’ 

(Ob  -  5) 

These  equations  have  been  verified  ty  experiments.  Since  the  calculation  procedure 
deaorlbed  in  this  chapter  is  baaed  on  more  simplifying  assumptions  than  the  results 
discussed  for  a  eonatant  wall  temperature*  it  seems  advantageous  to  determine  by  the 
method  of  this  chapter  only  a  correction  factor  which  has  to  be  applied  to  the 
constant  wall  temperature  heat  transfer  coefficients  in  order  to  take  into  account 
the  influence  of  a  wall  taaiperature  variation.  This  calculation  will  be  done  in 
the  following  lines: 

Neglecting  the  influence  of  the  variable  wall  temperature  on  heat  transfer  by 
using  the  heat  transfer  coefficient  for  a  constant  wall  temperature  the  heat  flow 
at  the  location  x  is  described  by  the  following  equation 


q  (x)  *  h(0,x)  [  T  (x)  -  T  1 
c  i  w  r  j 

Dram  this  and  the  previous  Btjwetlon  Cb  -  3  we  obtain  the  relation 


(Cb  •  6) 


,00  -  VVl>] 

h(0,x)  j^Tw(x)  -  TrJ 

(Cb  -  7) 
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for  the  ratio  of  the  actual  heat  flow  q(x)  to  the  heat  flow  q^x)  calculated  under 
the  ■implifylng  assumption.  for  a  ladur  boundary  l^rar  on  a  flat  plate,  thie 
equation  can  be  transformed  by  uee  of  tquntlon  Cb  -  4  into 


% (x)  TvM  ”  Tr  (Cb  -  8) 


and  for  a  turbulent  boundary  layer  on  a  flat  plate  Into  the  relationahip 

gw  m  Urr*  lTw(^>  - 

q0(x)  ^,(x)  -  Tr 


(Cb  -  9) 

Baeh  of  tbe  above  two  equation a  oan  be  aolved  ae  aoon  as  the  wall  tamparature  ae 

a  Auction  of  the  diatanoe  from  tbe  leading  edge  ia  preeoribed,  for  inetance  by 
numerical  or  graphical  integration  of  the  integrals  appearing  in  the  equations. 

Let  us  assume,  far  instanee,  that  the  wall  temperature  as  indicated  ia  the 
upper  diagram  of  Figure  13  is  prescribed.  This  temperature  variation  has  two 
discontinuities,  one  at  the  leading  edge  (  §•  ©  )  and  one  at  ft \  •  lbs  local 
heat  flow  at  the  location  x  for  laminar  flow  may  have  to  be  determined.  In  order 
to  find  the  Integral  ia  fquation  Cb  -  J  build  the  differential  quotient  dlW/df* 
multiply  this  value  with  and  plot  the  produce  over  §  ae  indicated 

in  the  lower  diagram  of  ngure  13.  The  area  under  this  curve  represents  the 
integral.  The  suamation  has  to  be  extended  over  O  and  f  (  an  3^  ; 

<^\<*V-VVi>]  •  [y°+>-Tr]  - 


(Ob  -  10) 

The  calculation  is  however  quite  tedious  because  the  integration  has  to  be  perfd— e< 
far  any  location  x  along  the  surface  for  which  the  local  heat  flow  is  to  be  determined. 
Additionally,  the  terns  under  tbe  integral  algae  tend  always  towards  infinity  when 
the  variable  &  under  the  lateral  sign  approaches  the  limit  x  as  shown  in  Figure  13. 
Therefore  the  contribution  of  the  region  near  x  to  tbs  integral  has  to  be  evaluated 
by  soon  special  procedure,  as  suggested  in  Appendix  XV* 

The  denoartnator  in  both  Iquatlons  Cb  -  8  and  Cb  -  9  oan  he  written  also  in  the 
following  formi 


(Cb  -  U) 


33 


It  is  in  tbs  example  equal  to  the  area  under  the  daahed  curve  in  Figure  13  plus  the 
sum  of  the  tvo  sudden  temperature  changes  at  the  leading  edge  and  at  sf),. 

The  procedure  described  assumed  that  the  vail  temperature  distribution  is 
known.  The  inverse  problem  to  find  the  wall  temperature  for  a  prescribed  oonveetive 
heat  flow  has  been  solved  in  Before nee  117.  Actually,  the  wall  temperature  has 
usually  to  be  determined  from  a  heat  balanoe  between  the  convective  heat  transfer 
and  the  radiative  heat  exchange  of  the  surface  with  the  surrounding  and  the  heat 
conducted  stray  into  the  solid  interior.  Generally,  this  problem  can  be  solved  only 
by  a  tedious  trial  and  error  procedure. 

In  this  situation  it  is  important  to  note  that  by  treatment  of  aoms  simplified 
case  it  has  been  determined  that  the  ratio  of  the  differences  between  wall  and 
stream  temperatures,  calculated  by  the  procedure  described  in  this  chapter) to  the 
corresponding  temperature  difference  based  on  a  simplified  procedure  which  neglects 
the  influence  of  the  wall  temperature  variation  on  the  heat  transfer  coefficient  is 
considerably  smaller  than  the  ratio  of  the  heat  transfer  coefficients  determined  by 
both  procedures. 

Li gh thill  (Ref.  119)  based  his  calculations  on  the  assumption  that  a  thin- 
walled  flat  plate  receives  heat  by  laminar  convection  from  a  gas  stream  and 
discharges  heat  by  radiation  into  a  surrounding  with  a  temperature  of  0°  Ranklne. 
Heat  conduction  within  the  wall  and  any  other  heat  removal  has  been  neglected. 

Using  the  expressions  given  by  equations  Cb  -  3  and  Cb  -  4  for  the  eonveotive  heat 
flow  into  the  surface  and  balancing  this  heat  flow  by  the  radiative  heat  transfer 
an  integral  equation  is  obtained  which  has  been  solved.  The  result  of  this 
calculation  which  is  contained  in  Figure  2  of  the  referenced  report  can  be  used  to 
make  a  comparison  with  a  calculation  based  on  heat  transfer  coefficients  far 
isothermal  surfaces.  It  has  been  found  that  the  wall  temperature  calculated  by 
Lighthill  differs  up  to  10  per  cent  from  the  absolute  wall  temperature  obtained  by 
the  simplified  calculations. 

Another  calculation  has  been  made  by  ftyson  and  Edwards  (Ref.  112)  for  the 
following  situation:  A  flat  plate  is  initially  at  a  constant  temperature  and  is 
suddenly  subjected  to  the  influence  of  a  high  velocity  gas  flow  of  constant  or 
exponentially  decreasing  velocity.  Under  the  influence  of  the  oonveetive  heat 
transfer  the  temperature  adjusts  itself  gradually  to  the  teaqierature  of  the  stream. 
Radiation  has  been  neglected  and  the  eonveotive  heat  flow  is  balanoed  with  the  heat 
stored  within  the  wall  (no  heat  conduction  within  the  wall  parallel  to  the  surface 
and  constant  temperature  throughout  the  wall  normal  to  the  surfaoe  are  assumed) 

The  use  at  Bquation  Cb  -  3  for  the  heat  flow  results  in  an  integro-diff ecenti si 
equation  which  has  been  solved  epd  compared  with  tbs  results  of  tbs  approximate 
prooedures  which  neglects  the  influence  of  e  temperature  variation  on  ths  hast 
transfer  coefficients.  It  was  found  that  ths  dlffsrenoe  between  wall  temperature 
and  static  stream  temperature  as  obtained  by  ths  sxaot  calculation  is  up  to  20  par 
oaat  largsr  than  tha  difference  obtained  by  the  Approximate  procedure.  A  difference 
of  10  per  cent  in  the  wall  temperature  as  found  1 y  Lighthill  or  at  20  par  oezxt  in 
ths  temper  store  difference  as  dstsrminsd  by  Biyson  corresponds  at  a  wall  tempera¬ 
ture  of  1500°  Rankins  and  a  atatle  stream  temperature  of  fcOO°R.  to  a  temperature 
difference  of  wound  200°F.  and  suoh  a  difference  mm  for  instance,  be  deciding  in 
the  decision  whether  a  certain  material  can  be  used  or  not  for  ths  Win  of  an 
aircraft.  It  m m  also  bs  that  ths  tsmpsreture  dlffsrenoe  on  special  locations  llfes 
ths  leading  sdgs  and  ths  trailing  sdgs  of  a  solid  wing  which  is  bains  basted  up  by 
aerodynamic  heating  (as  calculated  in  Rsfsrsnes  116)  mm  bs  largsr  than  in  the  above 
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examples.  Near  the  leading  edge  the  temperatures  will  be  increased  and  near  the 
trailing  edge  deereased  by  the  effeet  discussed  in  this  ohapter.  In  design 
calculations  it  will  therefore  hare  to  be  decided  from  the  situation  whether  the 
influence  of  a  variable  wall  temperature  has  to  be  included  into  the  calculations 
or  whether  the  much  simpler  procedure  which  neglects  this  influence  ia  considered 
as  sufficiently  accurate. 

A  check  on  the  ao curacy  of  the  method  discussed  in  thia  ohapter  has  been  made 
by  Lighthill  who  compared  it  with  Chapman  and  Rube  sin's  solution  for  laminar  flow 
and  a  wall  temperature  variation  according  to  Equation  Cb  -  1.  Agreement  within 
3  per  cent  was  obtained. 
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SECTION  V 


AUCTION  AND  HEAT  TRANSITS)  AT  LOW  DENSITIES 


The  solutions  for  friction  and  baat  transfer  which  hare  bean  presented  up  to 
now  are  subject  to  a  number  of  limitations.  One  of  these  limitations  will  be 
dismissed  in  the  present  chapter.  It  is  connected  with  the  fact  that  air  oouists 
of  single  molecules.  In  deriving  the  differential  equations  which  are  the  basis 
for  the  relationships  developed,  it  has  been  assumed  that  the  properties  of  the 
fluid  vary  steadily  throughout  the  field  investigated.  Actually  energy  and  momentum 
is  carried  without  change  ty  the  single  molecules  from  one  to  the  next  collision. 

As  long  as  the  mean  free  path  length  of  the  molecules  is  vary  small  compared  with 
the  length  dimensions  with  which  we  are  concerned,  we  can  assume  the  variation  to 
occur  steadily. 

The  mean  molecular  path  length  increases  inversely  proportional  to  the 
density  and  therefore  with  decreasing  density  it  has  to  be  expected  that  a  limit  is 
reached  where  the  free  molecular  path  length  is  of  the  same  order  of  magnitude  as 
the  characteristic  dimension.  The  ratio  of  free  molecular  path  length  to  the 
characteristic  length  is  called  Knudsen  number.  In  boundary  layer  flow  as  it  has 
been  considered  up  to  now  the  rate  of  change  of  the  properties  is  largest  in  a 
direction  normal  to  the  surface  throughout  a  length  equal  to  the  boundary  lqrer 
thickness,  for  this  type  of  flow  the  Knudsen  number  is  defined  as  the  ratio  c t 
mean  free  molecular  path  length  to  boundary  layer  thickness.  The  Knudsen  number 
Kh  may  be  expressed  by  the  flow  parameters  which  have  been  used  before,  namely  the 
Reynolds  number  and  the  Mach  number.  Aram  gas  kinetic  considerations  the  following 
relationship  is  derived 
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in  which  X !■  the  ratio  of  specific  heats  at  constant  pressure  and  volume. 

Ibr  boundary  layer  flow  the  Reynolds  number  is  based  on  the  boundary  layer 
thiekneaa  aa  characteristic  length.  Changing  to  the  Reynolds  number  based  on  the 
distance  along  the  surface  measured  from  the  leading  edge  results  in  the  following 
relationships  for  laminar  flow  along  a  flat  plate 


and  for  turbulent  flow 
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According  to  the  above  discuasion  the  molecular  structure  makes  itself  felt 
when  the  Khudsen  number  exceeds  a  certain  limit.  This  limit  can  be  well  illustrated 
in  a  diagram  which  uses  the  Reynolds  number  and  the  Mach  number  as  abscissa  and 
ordinata  respectively.  Such  a  diagram  is  presented  in  Figure  14 •  Inserted  in  the 
diagram  is  the  approximate  limit  between  turbulent  and  laminar  flow.  This  limit  is 
determined  by  a  large  number  of  factors  which  will  be  discussed  in  more  detail  in  a 
later  chapter.  The  line  given  in  Figure  14  should  be  taken  only  as  a  very  approxi¬ 
mate  indication.  The  line  with  the  break  on  this  boundary  between  laminar  and  tur¬ 
bulent  flow  indicates  the  conditions  for  which  the  Knudsen  number  for  laminar  flow 
as  an  expression  for  the  ratio  of  free  molecular  path  length  to  boundary  layer 
thickness  has  the  value  O.OO3.  In  the  field  on  the  right  hand  side  of  this  line  the 
free  molecular  path  length  is  less  than  1  per  cent  of  the  boundary  layer  thiokness 
and  the  relationships  developed  in  the  preceding  chapters  apply. 

With  increasing  distsnce  from  the  limiting  line  towards  the  left  hand  side  the 
free  molecular  path  length  increases  more  and  more  and  it  has  to  be  expeoted  that 
the  molecular  structure  of  the  air  causes  larger  and  larger  deviations  from  the 
conditions  encountered  in  a  continuum.  Near  the  limiting  line  these  deviations 
occur  only  within  the  immediate  neighborhood  of  the  solid  surface.  Their  sain 
effect  is  to  create  finite  flow  velocities  right  at  the  wall  surface  and  a  sudden 
variation  of  the  temperature  (temperature  Jump).  Accordingly  the  regime  Just  to 
the  left  of  the  limiting  line  is  called  slip  flow  regime. 

With  decreasing  density  finally  a  state  is  reached  where  now  inversely  the 
free  molecular  path  length  is  such  larger  than  any  dimension  of  the  body*  In  such 
a  case  the  molecules  impinge  on  the  body  surface  practically  with  the  momentum  and 
energy  which  prevail  in  large  distance  from  the  body.  The  field  in  Figure  14  in 
which  this  condition  applies  is  termed  free  molecular  flow  regime.  The  limit  for 
this  regime  is  usually  fixed  by  the  condition  that  the  value  M^Re  which  is  not 
much  different  from  the  Knudsen  number  has  the  value  10  (Ref.  140) 


Ma  /  Re  -  10 
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No  boundary  layer  exists  under  this  condition.  Therefore  Re  is  based  on  a  character¬ 
istic  body  dimension.  The  line  described  by  Equation  D  -  4  1*  also  shown  on  the 
left  hand  side  of  Figure  14 . 


Another  limit  to  the  relationship  for  contininuum  flow  presented  in  the 
previous  chapters  is  connected  with  the  fact  that  with  decreasing  Reynolds  number 
and  increasing  Mach  number  the  boundary  layer  thickness  increases.  The 
layer  equations  on  which  the  solutions  are  based  become  invalid  when  the  boundary 
layer  thiokness  exceeds  a  certain  limit.  The  displaoeamnt  of  the  flow  by  the 
boundary  layer  becomes  then  so  large  that  it  onuses  the  pressure  to  vary  notably 
along  a  flat  plate.  This  effect  is  especially  large  where  the  flow  field  is  limited 
by  a  shock  ware  and  when  this  shook  wave  is  comparatively  near  to  tbs  surface. 
Different  authors  have  made  estimates  on  the  limit  beyond  which  this  boundary  layer 
shock  wave  interaction  becomes  important.  The  results  indioate  generally  that  the 
limit  is  situated  near  to  the  limit  between  the  slip  flow  and  continuum  flow  region 
soswwhat  displaced  toward  the  left  for  low  Mach  numbers  and  towwd  the  right  at  high 
Mach  numbers.  Correspondingly,  for  high  Itoeh  numbers  the  slip  effect  will  be  the 
first  one  which  modifies  the  continuum  flow  conditions  whsn  the  Reynolds  number 
decreases,  whereas  at  low  Mach  numbers  the  mentioned  boundary  layer  shook  wave 
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interaction  will  cause  the  first  modification. 


In  general  the  friction  and  heat  transfer  conditions  outside  the  range  of 
continuum  flow  are  extremely  complicated  and  therefore  little  investigated  analyti¬ 
cally  as  well  as  experimentally.  A  comparatively  good  knowledge  exists  only  for  the 
regime  far  to  the  left  in  Figure  14  which  has  been  termed  'free  molecular  flow* 
regime.  Very  such  work  remains  to  be  done  in  the  'slip  flow*  regime  and  especially 
in  the  'mixed  flow*  regime. 

In  this  situation  it  is  important  to  visualize  the  regions  in  the  Reynolds - 
Mach  number  field  which  will  be  of  special  importance  in  future  aeronautical 
applications.  In  this  paper  we  are  concerned  with  cooling  problems  and  from  this 
point  of  view  we  can  obtain  the  desired  information  in  the  following  ways  Cooling 
problems  will  always  become  critical  when  the  skin  of  the  aircraft  exoeeds  certain 
temperature  limits.  The  surface  temperature  of  the  airoraft  is  determined  by  the 
combined  effects  of  eonveetive  heat  transfer  from  the  atmosphere  to  the  surface,  of 
emission  of  radiation  from  the  surface,  of  irradiation  from  the  sun,  and  finally  of 
cooling  either  by  storage  of  heat  in  the  solid  material  of  the  airoraft  or  by  soma 
cooling  process.  All  of  these  heat  exchange  processes  depend  on  the  flight  altitude 
and  the  flight  velocity,  both  of  which  can  be  expressed  by  the  Reynolds  number  and 
Mach  number  at  which  the  flight  occurs .  Two  curves  calculated  on  the  basis  of  the 
mentioned  considerations  in  Appendix  II  have  been  inserted  in  Figure  14 .  Along 
these  curves  a  wall  temperature  of  1500°R  is  obtained  by  the  surface  of  the  aircraft 
when  the  surface  radiates  heat  like  a  black  body  and  when  cooling  by  heat  storage  or 
by  some  cooling  process  does  not  occur.  One  of  the  curves  holds  for  a  looation  on 
the  aircraft  which  haa  a  distance  of  one  foot  from  the  leading  edge  and  the  second 
curve  for  a  distance  of  5  ft.  As  break  in  the  curve  at  Re  ■  2-10®  is  connected  with 
the  fact  that  cooling  conditions  are  more  favorable  in  a  laminar  boundary  layer  than 
in  a  turbulent  one.  Depending  on  the  specific  conditions  the  looation  of  the 
transition  point  may  shift  to  the  right  or  left  from  the  position  indicated  in  the 
figure. 

Reynolds  number  and  Mach  number  together  with  the  other  data  which  have  been 
mentioned,  fix  the  flight  altitude  along  the  two  curves  and  the  corresponding  values 
are  also  indicated.  The  left  hand  curve  holds  also  for  tbs  condition  that  tbs 
distance  L  is  equal  to  5  ft.  but  the  emislvity  of  the  surface  for  heat  radiation  la 
equal  to  0.2.  The  two  curves  therefore  indicate  the  Reynold  and  Maoh 
conditions  far  which  an  uncooled  surface  of  an  airoraft  assumes  a  temperature  of 
1500°R  in  steady  flight.  If  cooling  is  applied  to  the  surface  either  by  heat 
storage  which  is  possible  for  short  flight  duration  or  by  some  oooling  process  then 
the  flight  velocity  and  correspondingly  the  Mach  number  can  be  higher  fer  a  certain 
Reynolds  number  without  increasing  the  surface  temperature  above  the  value  of 
1500°R.  For  this  condition  therefore  the  flight  would  occur  at  a  point  which  Is 
located  in  Figure  14  above  and  to  tbs  right  of  ths  curve.  Aircraft  or  mi  sails 
during  their  flight  m ay  encounter  conditions  as  given  by  their  Reynolds  and  Mech 
number  which  belong  to  ths  range  located  below  and  to  tbs  left  of  ths  two  curves  in 
Figure  14.  This  then  indicates  that  ths  wall  temperature  even  for  steady  flight 
will  be  lower  than  the  assumed  limit  of  1500°R  end  generally  that  euoh  flight  condi¬ 
tions  ora  not  critical  from  a  cooling  standpoint.  Ths  conclusion  therefore  is  that 
the  range  near  the  two  curves  in  Figure  14  is  ths  one  for  which  infonaatlon  on  hast 
transfer  is  most  important.  With  this  in  view  we  will  new  review  ths  information 
on  friction  and  heat  transfer  which  is  available  for  ths  regimes  outside  ths 
continuum  flew  region. 
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Comparatively  good  information  is  available  for  the  conditions  in  the  free 
molecular  flow  regime.  The  kinetic  theory  of  gases  gives  the  possibility  to 
calculate  friction  and  heat  transfer  to  objects  in  a  flow  with  uniform  velocity  and 
uniform  temperature,  figure  15  and  16  show  such  data.  The  figures  have  been  taken 
from  Reference  135.  Figure  15  gives  Stanton  numbers  for  a  monatomic  and  a  diatomic 
gas.  The  information  on  a  monatomic  gas  is  important  because  the  air  molecules 
dissociate  at  very  high  altitudes.  The  value  0(  by  which  the  Stanton  numbers  are 
divided  in  the  figure  is  the  accomodation  coefficient.  Its  value  for  air  is  not 
too  well  known.  Usually  it  is  assumed  to  be  equal  to  0.9.  An  interesting  fact  may 
be  observed  in  Figure  16  which  presents  the  recovery  factors  for  two-molecular  and 
one-molecular  gases  as  a  function  of  the  Mach  number.  It  can  be  seen  that  quite  a 
number  of  configurations  among  them  the  flat  plate  in  a  flow  parallel  to  its 
surface  have  recovery  factors  which  are  larger  than  one.  We  have  therefore  to 
expeot  that  with  decreasing  density  the  recovery  factor  which  in  continuum  flow  of 
air  has  a  value  of  O.85  or  0.9  increases  and  assumes  values  above  1  in  the  free 
molecular  range.  This  trend  has  actually  been  verified  experimentally. 

Information  on  heat  transfer  in  the  mixed  flew  regime  is  very  limited.  This 
is  an  unfortunate  fact  sinoe  we  have  seen  that  this  region  is  of  special  importance 
for  future  aeronautical  applications.  Calculations  which  have  been  made  rest  on  a 
comparatively  unsecure  basis  sinoe  it  is  today  even  not  known  exactly  what  the 
differential  equations  are  which  describe  the  flow  and  energy  transfer  process  in 
this  region.  Experimental  information  has  been  obtained  at  the  Asms  Laboratory  of 
the  IttCA  and  at  the  University  of  California  in  Berkeley.  Such  information  is  very 
valuable  in  view  of  the  difficulties  of  obtaining  analytio  information.  However, 
these  tests  are  restricted  to  Mach  numbers  below  4  and  from  Figure  14  it  can  be 
seen  that  the  vsry  high  Mush  number  condition  is  especially  Important  in  this  range. 
It  is  not  at  all  sure  whether  information  obtained  at  small  Mach  numbers  can  be  used 
for  the  high  Mach  number  regime  as  well.  One  additional  factor  a my  be  observed  in 
Figure  14.  Any  information  obtained  from  analysis  and  all  the  experimental 
information  obtained  so  far  deal  with  laminar  flow  region.  The  curves  indicating 
the  flight  conditions  under  which  the  surface  of  an  aircraft  assumes  a  temperature 
of  1500°R  however.  Indicate  that  at  least  part  of  the  flight  of  aircraft  mid 
missiles  has  to  bs  expected  to  be  in  the  turbulent  region,  at  least  near  the  slip 
flow  regime  and  there  exists  practically  no  knowledge  on  heat  transfer  under 
turbulent  flow  conditions  outside  of  the  oontinuum  regime. 

In  the  rest  of  this  chapter  information  will  be  discussed  whioh  has  resulted 
from  investigations  in  the  slip  flow  and  the  mixed  flow  regime  for  1  — flow 
conditions.  Hear  the  oontinuum  regime  the  normal  boundary  layer  equations  still 
hold  to  a  good  approximation.  Only  the  boundary  conditions  have  to  be  changed  to 
take  into  account  the  finite  velocities  and  the  sudden  change  in  temperature 
00 curing  on  e  solid  surface  (slip  flow  and  temperature  Jump)  and  the  interaction 
between  shock  waves  and  boundary  layer.  Meslen  (Ref.  132)  has  mads  a  calculation 
on  thla  baaia  whioh  is  essentially  a  perturbation  of  tbs  normal  oontinuum  solution. 
This  calculation  determines  to  a  first  approximation  the  changes  which  occur  in  tbs 
norad.  oontinuum  solutions  when  the  flew  enters  the  slip  flow  regime.  These  results 
of  the  calculations  are  expressed  by  the  following  equations  for  tbs  skin  friction 
coefficient 
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and  for  the  heat  transfer  coefficient 
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In  both  equations  the  subscript  zero  indicates  the  normal  continuum  solution,  for 
instance,  as  discussed  in  the  previous  chapters  of  this  paper.  The  subscript  1 
indicates  the  change  which  this  continuum  solution  undergoes  to  a  first  approximation 
and  the  parameter  £  is  essentially  the  Khudaen  number: 
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The  calculation  procedure  is  restricted  to  tne  assumption  of  a  constant  specific 
heat  for  the  air  but  is  applicable  to  any  variation  of  visoosity  and  irandtl  number 
with  temperature .  Tsble  II  gives  the  results  of  a  numerical  evaluation.  The 
friotlon  factors  and  heat  transfer  coefficients  as  calculated  by  Young  and  Janssen 
(Ref.  59)  are  used  as  the  continuum  solution  and  the  changes  oc curing  to  a  first 
approximation  have  been  calculated  by  the  procedure  described  in  Reference  132.  from 
lible  II  and  Aram  the  Kquatlons  D-5toD-7it  may  be  observed  that  the  shearing 
strees  inereasea  above  the  value  in  continuum  flow  when  the  slip  flow  region  is 
entered.  This  Increase  is  caused  by  the  boundary  layer  shock  wave  interaction.  The 
heat  transfer  paranmter  decreases  below  the  value  in  oontinuum  flow  as  a  consequence 
of  the  temperature  Jump  00 curing  at  the  surface  under  slip  flaw  conditions.  Both 
changes  are  of  the  order  of  magnitude  of  1  per  cent  at  the  limit  which  is  indicated 
in  Figure  14  between  continuum  flow  and  slip  flow  regime.  The  difference  between 
the  actual  and  the  continuum  solution  increases  with  decreasing  Hnudsen  number  and 
is  about  10  per  cent  when  the  par  master  Me/  assumes  the  value  0.1  indicating 
that  the  Area  molecular  path  length  has  increased  to  approximately  1/10  of  the 
boundary  layer  thickness. 

An  indication  of  how  friction  and  heat  transfer  change  throughout  the  whole 
region  between  oontinuum  and  free  molecular  flow  any  be  obtained  from  Figures  17 
and  18  which  are  taken  Aram  a  paper  by  Shaaf  (Ref.  I38)  susmarising  the  work  of 
the  Berkeley  group.  The  figures  contain  as  dashed  lines  the  solutions  obtained 
Aram  the  boundmry  layer  equations  and  valid  in  the  oontinuum  flow  regime  as  well  as 
the  solutions  obtained  for  the  free  molecular  flow  regime.  The  full  line  gives  the 
reoaaamnded  correlation.  Some  experimental  points  have  also  been  included.  The 
Reynolds  and  Meoh  numbers  appearing  in  Figure  18  are  based  on  the  conditions  behind 
the  shook  wave.  The  value  #  is  half  the  apex  angle  of  the  cone. 
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SECTION  VI 


FRICTION  AND  HEAT  TRANSFER  UNDER  HICK  TEMPERATURE  CONDITIONS 


Id  the  flight  of  aircraft  through  the  atmosphere  with  extreme  Maoh  number  very 
high  temperatures  are  encountered  within  the  boundary  layer.  These  temperatures 
cause  the  air  to  dissociate  and  finally  to  ionize  (at  temperature  above  20,000°R). 

In  this  chapter  the  limitations  will  be  discussed  which  are  imposed  on  the  boundary 
layer  solutions  in  the  previous  chapters  under  high  temperature  conditions. 

Figure  19  has  been  prepared  to  indicate  what  temperatures  have  to  be  expected 
within  the  boundary  layer  at  high  Maoh  numbers.  The  figure  holds  far  a  stream 
static  temperature  of  400°R*  When  the  surface  assumes  the  recovery  temperature 
then  this  temperature  will  also  be  the  highest  static  temperature  rhich  is  encount¬ 
ered  within  the  boundary  layer.  This  recovery  temperature  of  a  flat  plate  ia 
indicated  in  the  figure  for  laminar  and  turbulent  flew  conditions.  The  full  line 
has  been  calculated  assuming  constant  specific  heat  and  a  recovery  fectar  of  0.85 
for  laminar  flow.  The  dashed  line  is  calculated  for  a  recovery  faotor  of  0.9 
corresponding  to  turbulent  flow  and  again  assuming  constant  specific  heat.  A 
third  curve  has  been  added  from  the  results  of  Reference  134,  which  considered  the 
specific  heat  to  vary,  in  order  to  indicate  what  variation  in  the  recovery  tempera¬ 
ture  is  caused  by  the  variation  of  the  apeoific  heat. 

In  aircraft  applications  the  surface  temperature  will  usually  be  cooled  by 
radiation  or  by  some  cooling  process.  Under  such  conditions  the  maximum  temperature 
occurs  somewhere  within  the  boundary  layer  (See  also  Fig.  5)*  Curves  have  been 
added  in  the  figure  which  indicate  this  maximum  temperature  within  the  boundary 
layer  for  the  case  that  the  wall  temperature  is  cooled  down  to  the  streets  static 
temperature  and  for  the  case  that  the  absolute  wall  temperature  is  four  times  the 
free  static  temperature  (at  a  temperature  of  l600°R).  Jbr  a  gas  with  a  frandtl 
number  equal  to  1  it  can  be  easily  derived  thet  the  difference  between  the 
temperature  within  the  bovmdary  layer  and  the  stream  static  temperature  is  equal 
to  one  fourth  of  the  difference  between  recovery  temperature  and  statio  temperature 
in  the  stream  for  laminar  as  well  as  far  turbulent  flow  under  the  condition  that  the 
wall  temperature  is  equal  to  the  static  stream  temperature.  For  the  calculation  of 
the  curves  in  Figure  19  it  has  been  assumed  that  this  condition  still  holds 
approximately  far  air  with  a  Brandt 1  number  of  0.72. 

It  can  be  determined  from  Reference  165  that  dissociation  starts  at  an  altitude 
of  5°  miles  in  the  air  with  a  temperature  around  3,000°R.  and  at  sea  level  with  a 
temperature  around  6,000°R.  The  Figure  19  in  which  these  limits  are  also  inserted 
indicates  that  no  dissociation  within  the  boundary  layer  is  expeoted  up  to  Maoh 
numbers  of  6  or  10  for  the  case  where  the  wall  temperature  is  kept  at  the  recovery 
temperature  and  that  dissociation  will  occur  only  at  Mach  numbers  above  12  to  18 
when  the  wall  temperature  is  cooled  to  values  below  16009R.  The  latter  case  is  the 
one  which  is  more  likely  to  occur  on  aircreft. 

Two  papers  (Ref.  142,  and  143)  have  been  published  recently  whioh  investigate 
what  effect  dissociation  will  have  on  friction  and  heat  transfer  in  laminar  flow 
along  a  flat  plate,  Such  calculations  are  connected  with  considerable  uncertainties. 
They  are  for  instance  based  on  the  assumption  that  the  air  adjusts  Itself  at  each 
location  within  the  boundary  layer  to  the  equilibrium  dissociation  belonging  to  the 
local  temperature.  It  is  actually  quite  doubtful  whether  the  fluid  particles  stqr 


for  a  sufficient  time  in  the  high  temperature  zone  to  assume  equilibrium  dissocia¬ 
tion.  Also  diffusion  of  the  different  gases  into  which  the  air  dissociates  into 
the  cooler  regions  may  occur.  All  these  facts  are  neglected  in  the  calculations. 
Another  uncertainty  is  based  on  our  insufficient  knowledge  of  property  values  of  a 
gas  at  extremely  high  temperature  and  under  the  condition  of  dissociation.  In  this 
point  the  assumptions  made  in  both  sets  of  calculations  differ  considerably,  for 
instance  Moore  (Ref.  143)  assumes  that  the  fraadtl  number  in  diasooiated  air 
increases  up  to  values  around  6  whereas  Crown  (Ref.  142)  reasons  that  even  with 
dissociation  the  frandtl  number  can  never  assume  values  larger  than  1.  Nevertheless 
both  calculations  cone  in  some  important  points  to  the  same  conclusion.  It  ia  found 
that  dissociation  under  the  assumed  conditions  has  a  negligible  effect  on  the  skin 
friction  and  on  the  heat  flow  entering  or  leaving  the  surface  per  unit  area  as 
long  as  the  s  urface  temperature  is  kept  at  values  which  are  below  the  temperature  at 
which  the  dissociation  in  the  air  starts.  Ibis  situation  holds  even  though  very 
much  higher  temperatures  may  be  encountered  within  the  boundary  layer.  The  condition 
that  the  wall  temperature  is  lower  than  the  dissociation  temperature— lower  than 
say  3.000°R— will  be  fulfilled  in  almost  all  aeronautical  applications. 

Crown  finds  additionally  in  his  calculations  which  are  based  on  the  assumption 
that  the  frandtl  number  is  constant  and  that  the  product  of  viscosity  and  density 
is  a  linear  function  of  the  enthalpy  that  the  recovery  factor  can  be  represented 
within  a  few  per  cent  by  the  well-known  relationship 


r±  -  YR 

(E  -  1) 

He  also  finds  in  Kquation  38  of  his  report  an  expression  which  can  be  transformed 
into  the  following  relation  fur  the  ratio  of  Stanton  numbers  to  friction  factor 


®f  Tr  -  T« 

(E  -  2) 

This  equation  is  identical  with  Equation  fla  -  8  and  indicates  that  the  Stanton 
number  is  practically  Independent  of  dissociation  when  baaed  on  enthalpy. 

Moore  performed  his  calculation  considering  that  the  frandtl  number  for  air 
varies  with  temperature.  Enthalpy  recovery  factors  calculated  from  the  recovery 
temperatures  contained  in  his  report  have  been  compared  with  enthalpy  recovery 
factors  calculated  from  equation 


•  yur* 


(*  -  3) 


with  the  frandtl  number  introduced  at  the  reference  enthalpy  as  given  by  Iqiation 
Be  -  26.  The  agreement  is  not  too  good  for  this  set  of  data,  formula  (E  -  3) 
gives  recovery  temperatures  ehich  are  up  to  20  per  cent  too  high,  for  design 
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calculation*  the  actual  recovery  temperature  under  conditions  of  dissociation  is  of 
little  interest  since  it  is  so  high  that  no  presently  known  material  can  withstand 
it.  Surfaces  of  aircraft  flying  at  such  Mach  numbers  will  hare  to  be  cooled 
considerably.  Into  calculations  of  the  oonreotire  heat  transfer  under  such 
conditions  the  recovery  temperature  enters  only  as  a  means  to  express  the  heat  flow 
conveniently  (Squat ion  Ba  -  3)  and  for  this  purpose  the  reoorery  temperature 
determined  by  Equation  £  -  3  gives  better  results  then  the  actual  recovery 
temperature. 

Crown  presented  in  his  paper  formulas  for  a  calculation  of  friction  and  heat 
transfer  which  are  valid  for  any  temperature  dependence  of  the  properties  in  the 
same  w ay  as  the  equations  presented  in  the  proceeding  chapters.  The  method 
discussed  in  this  paper  is  believed  to  have  the  advantage  that  it  is  in  accord  with 
the  procedures  with  which  in  heat  transfer  calculations  a  variation  of  properties 
is  accounted  for. 

The  findings  which  have  been  discussed  are  very  useful  information  for  an 
estimate  of  heat  transfer  coefficients  under  high  temperature  conditions.  Of 
course,  it  has  to  be  kept  in  mind  that  they  are  based  on  assumptions  which  still 
have  to  be  booked  up  by  experiments  and  that  they  hold  directly  only  for  laminar 
boundary  layers.  Experimental  investigations  under  high  temperature  conditions  are 
needed  for  more  accurate  knowledge  of  friction  and  heat  transfer  conditions,  ftiey 
are  however,  extremely  difficult  to  perform.  In  wind  tunnels  for  Instance,  it  is 
almost  hopeless  to  increase  the  temperatures  to  such  values  that  on  surfaces,  which 
are  cooled  to  conditions  similar  to  the  ones  expected  on  high  velocity  aircraft  the 
influence  of  dissociation  can  be  studied,  from  Figure  19  for  example,  it  can 
easily  be  determined  that  for  a  static  temperature  of  400^.  in  an  air  stream  and  a 
pressure  corresponding  to  an  altitude  of  30  miles  the  total  temperature  or  stagnation 
temperature  has  to  be  raised  to  values  above  8,000  to  10,000°R.  before  aay 
dissociation  will  occur  in  boundary  layers  which  develop  on  surfaces  cooled  to  a 
temperature  of  l600°R.  It  is  to  be  expected  that  even  considerably  higher  total 
temperatures  are  necessary  to  obtain  a  marked  effeot  of  dissociation  on  frietion 
and  heat  transfer  at  the  surface. 
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SECTION  m 


TRANSITION  IRON  LAMINAR  TO  TURBULENT  HOW 


In  order  to  apply  the  information  presented  in  the  previous  chapters  it  is 
necessary  to  have  additional  knowledge  as  to  which  portion  of  a  surface  is  oovered 
ty  turbulent  boundary  layers  or,  in  other  words  where  the  transition  from  laminar 
to  turbulent  flow  occurs.  This  point  is  the  weakest  link  in  the  whole  procedure 
to  predict  wall  temperatures  and  heat  transfer  coefficients.  The  present  day 
knowledge  has  been  recently  sumnarised  by  Qezley  in  Reference  14?.  The  reader 
therefore  is  referred  to  this  reference.  Only  one  figure  may  be  added  here  which 
presents  the  results  of  a  stability  calculation  by  Van  Erlest  (Ref.  145)  sinoe  it 
extends  to  considerably  higher  Mach  numbers  than  the  information  contained  in 
Reference  147*  figure  20  shows  as  a  funetion  of  Mach  number  and  of  wall  to  stream 
atatic  temperature  the  Reynolds  number  at  which  boundary  layera  on  a  flat  plate 
bacons  unstable  for  small  fluctuations  of  certain  ware  length. 

The  full  lines  in  the  figure  indicate  the  critical  Reynolds  numbers  for  air  of 
394?R.  static  stream  temperature  assuming  a  JTandtl  number  0.72  and  a  viscosity 
which  follows  Sutherland's  relation.  Limited  by  the  curve  Re  ■  ee  is  a  field  in 
which  the  flow  fluctuations  always  damp  out.  The  border  of  this  field  wee  aleo 
calculated  for  a  gas  with  a  viscosity  varying  proportional  to  temperature  (  fM-  * 
constant)  and  for  two  IVandtl  numbers.  The  dashed  line  holds  far  these  conditions. 
It  stay  be  recognised  that  both  parameters  have  considerable  influence  on  stability. 

It  is  known  especially  from  experiments  at  a  lower  Mach  number  that  a  certain 
time  or  length  is  necessary  far  these  fluctuations  to  increase  to  such  an  amount 
that  they  transform  the  flow  to  turbulence.  Accordingly  the  critical  Reynolds 
number  for  transition  is  by  a  factor  up  to  100  larger  than  the  Reynolds  number 
predicted  from  figure  20.  However,  qualitatively  the  trends  predicted  by  stability 
oaloulations  have  been  verified  generally  by  experiments.  The  figure  can  therefore 
be  used  to  make  estimates  on  the  influence  of  Mech  number  and  wall  to  stream 
temperature  ratio  on  the  critical  Reynolds  number  at  which  transition  occurs.  Iren 
in  this  respect  however,  one  hae  to  be  careful,  sinoe  Oesely  presents  experimental 
information  which  in  some  oases  does  not  conform  with  predictions  of  the  stability 
theory,  figure  20  together  with  the  information  in  Reference  147  has  been  used  to 
estimate  the  orltioal  Reynolds  number  for  a  wall  to  stream  temper,  ture  ratio  equal 
to  3  and  this  estimated  curve  has  been  entered  into  Figure  14*  W  oompering  this 
curve  with  the  two  lines  whieh  describe  the  flight  conditions  for  on  aircraft  with 
a  surface  temperature  of  1500°R.  it  mqr  be  concluded  that  probably  up  to  an  altitude 
around  30  or  40  miles  the  boundary  layers  will  be  turbulent  over  the  major  portion 
of  the  surface.  At  higher  altitudes,  when  the  alroraft  enters  the  mixed  flow 
region  and  moves  towards  the  free  moleoular  flow  regime,  the  turbulence  has  to  die 
down  somewhere  because  certainly  no  turbulence  oan  be  imagined  in  the  free  moleoular 
flow  regime. 


AlttKWnOM 


£.  ft.  Van  Driest  recently  extended  hie  calculations  of  friction  and  beat 
transfer  in  a  laminar  boundary  layer  on  a  flat  plate  to  oover  a  rery  wide  range  of 
Mach  numbers  (Ma  0  to  16)  of  total  stream  temperatures  (T*  0  to  2000°/)  and  of  the 
ratios  enthalpy  at  surface  to  enthalpy  in  stream  (i^/i  X  to  6) (ftef .  182) 
Additionally,  results  for  400°R  stream  temperature  are  presented. 

lhe  method  to  calculate  friction  factors,  recovery  factors  and  heat  transfer 
paraasters  2  St/c-  with  the  reference  temperature  proposed  in  this  report  has  been 
checked  against  the  results  of  Van  Driest 'a  calculations.  Agreement  of  the 
friction  factors  was  within  4  per  cent,  of  the  enthalpy  recovery  factors  within 
2  1/2  per  cent  and  of  the  heat  transfer  parameter  within  1  per  oent. 
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APPENDIX  I 


REFERENCE  TEMPERATURE 


In  this  appendix  it  will  be  described  how  the  Equation  Ba  -  16  has  been 
obtained  from  the  data  in  Reference  59  •  In  this  paper  the  reference  temperature  T' 
is  tabulated.  This  temperature  is  determined  in  such  a  way  that  the  friction 
factor,  calculated  from  the  boundary  layer  solutions  and  based  on  property  values 
at  the  Temperature  T'  agrees  with  the  friction  factor  far  a  constant  property  value 
field. 

The  temperature  T*  will  now  be  approximated  by  a  temperature  T*  as  gives  by  an 
equation  of  the  form: 


T»  -  ?,♦  A< V  Tg)  ♦  B(Tr-  T#) 

(I  -  1) 

Tables  2,  3,  and  4  in  Reference  59  represent  the  temperature  T'  for  the  condition 
that  the  wall  temperature  is  equal  to  the  recovery  temperature.  In  this  oase 
Equation  1*1  changes  to: 

T*A,  -  1  ♦  (A  ♦  B)  (V*#  -  1) 

(I  -  2) 

This  equation  is  represented  by  a  straight  line  in  Figure  21.  In  this  diagram  the 
temperature  ratio  T'/T.  taken  from  Reference  59  is  plotted  over  Ty/Tg  .  It  can  be 
obeerved  that  the  points  representing  the  exact  values  of  the  reference  temperature 
can  be  well  approximated  by  a  straight  line.  The  tangent  of  the  line  yields: 

A  ♦  B  »  0.72. 

Tables  5.  6,  7  in  Reference  59  contain  the  reference  temperature  T'  for 
conditions  where  the  wall  temperature  is  different  from  the  recovery  temperature. 

In  order  to  obtain  the  constants  A  and  B  individually  from  these  data.  Equation 
I  -  1  is  brought  into  the  fora: 


T*  -  T.  -  A(TW  -  T,)  ♦  (o.72  -  A)  (Tr  -  T,)  (j  _  3) 
or: 


0,72  fl  ■  £  -  i ;  > 

(I  -  4) 

This  relation  is  again  represented  by  a  straight  lias  in  Figure  22.  In  this 
diagram  are  plotted  two  values  0.72  (T  -T»)/Tg  over  (T-T  )/T  .  These  data  were 
calculated  from  the  data  in  Reference  59  •  It  is  seen  that  these  values  can  also  be 
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wall  approximated  by  a  straight  11m.  Dm  tangent  of  the  inserted  11m  gives 
A  =  0.50.  which  yields:  B  =  0.22.  In  this  way  Iquation  Ba  -  16  is  established. 

The  sane  procedure  was  carried  out  replacing  in  all  of  the  above  equations 
the  temperature  by  the  enthalpy  corresponding  to  the  respective  temperature.  The 
points  indicated  by  the  white  symbols  in  Figures  21  and  22  were  obtained  and  it  may 
be  obeerved  that  the  seam  straight  11m  approximates  the  enthalpy  ratios  as  well  as 
the  temperature  ratios. 
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APPENDIX  II 


HEAT  BALANCE  EQUATION 

The  temperature  which  the  skin  of  an  aircraft  assumes  is  to  be  calculated  for 
the  following  conditions!  The  aircraft  flies  horizontally  with  constant  speed  for 
a  time  sufficient  for  the  skin  to  assume  its  steady  state  temperature.  No  cooling 
process  is  used  to  remove  heat  from  the  skin  and  the  temperature  is  determined  by  a 
balance  between  the  heat  exchanged  with  the  universe  by  radiation.  The  skin  surface 
may  be  considered  "black"  with  respect  to  radiation. 

The  heat  balance  per  unit  area  is  written* 

*4.  (*r  “  “  ^eun 


(II  -  1) 

For  turbulent  flow  the  Equations  Ba-38,  Ba-7,  and  Ba-8  are  used  to  express  the 
heat  transfer  coefficient  by  the  Reynolds  number.  The  recovery  enthalpy  is 
expressed  by  Equation  Ba-21,  Ba-47  and  the  relation* 


-  t=W  i. 

(II  -  2) 

this  results  in* 

0.0296  («**)0,8  (ijl  ♦  0.18  (Ma)2]  - 

(II  -  3) 

This  equation  gives  for  the  prescribed  wall  temperature  of  Tw  =  1500°R  and  stream 
temperature  the  necessary  relation  between  Reynolds  -  and  Mach  numbers.  The  stream 
temoerature,  however,  is  connected  with  the  flight  altitude  and  this  value  is  a 
'‘unction  of  Re  u«d  to  through  th~  following  relation  for  the  density  of  airs 


_JL  Si 

x  a,  Ma 


(II  -  4) 


In  addition,  Pr  and  fx  d-vv^nd  on  the  reference  ,»nihalpy  and  therefore  on  stream 
temoerature  (temperuture  of  atmosphere  in  altitude)  and  on  the  Mach  number.  The 
problem  to  determine  the  Reynolds  number  which  corresponds  to  a  prescribed  Mach 
number  may  therefore  be  solved  by  a  trial  and  error  procedure*  Assune  a  stream 
temperature,  calculate  from  E  iation  (II  -3)  the  Reynolds  number,  from  Equation 
(II  -  4)  the  air  density  and  jheck  in  tables  for  properties  of  the  atmosphere 
whether  the  assvaed  stream  temperature  corresponds  to  the  altitude  as  determined  by 
the  air  density.  The  result  of  such  calculation  for  x=l  and  x»5  ft  is  inserted  in 
Figure  14. 
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Laminar  flow  conditions  occur  mainly  in  the  mixed  flow  regime  and  the  heat 
transfer  coefficients  for  this  region  would  have  to  be  taken  from  Figure  18.  This, 
however,  makes  the  calculation  v»ry  inconvenient.  Therefore  the  calculations  were 
made  for  laminar  continuum  flow,  expressing  the  heat  transfer  coefficient  in 
Equation  (II  -  1)  by  Equations  (Ba-37)  and  (Ba-28)  and  were  extended  only  as  far 
as  the  continuum  relations  are  expected  to  hold  still  with  a  reasonable  accuracy. 

The  curves  of  constant  wall  temperature  in  the  upper  right  hand  corner  were 
calculated  with  the  heat  transfer  coefficient  for  free  molecular  flow. 


49 


APPENDIX  IU 


LIMIT  TOR  ACCELERATION  FOR  WHICH  BOUNDARY  LATE R- 
nn*  CAN  EE  REGARDED  AS  QUASI -STEAIY 


F.  1-bore  shows  in  Ref.  5  that  boundary  layer  flow  along  a  flat  plate  may  be 
treated  as  quasi-steady  when  the  dimsnsionless  parameters: 


x  v;  x2  v«  x°  v<n) 
v'i  *  73  '  vn*I 

•  1  'i 


(hi  -  1) 


are  small.  In  these  parameters  V*  is  the  first  time  derivative  of  the  velocity. 

V*  the  second  etc.  Die  influence  of  the  first  one  of  these  parameters  on  the 
friction  factor  and  recovery  factor  of  an  insulated  flat  plate  with  laminar  boundary 
layer  has  been  evaluated.  It  was  found  that  the  shearing  stress  is  larger  by 
2.5  per  cent  and  the  recovery  factor  smaller  by  2.8  per  cent  than  the  steady  state 
values  for  flow  with  constant  acceleration  when  the  condition: 


(HI  -  2) 

is  fulfilled.  This  means,  for  instance,  that  the  acceleration  of  a  plate  moving 
with  1000  ft  per  sec.  speed  has  to  be  larger  than  10,000  ft/sec2  or  larger  than 
300  times  the  gravitational  acceleration  before  deviations  larger  than  2 .5  per  cent 
from  steady  state  conditions  are  encountered  at  a  location  of  1  ft  distance  tram 
the  leading  edge.  For  a  distance  of  10  ft  the  necessary  acceleration  is  1000 
ft/ sec2  or  30  g. 
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aKENDU  IV 


EVALUATION  OF  INTEGRALS  (Cb  -  8)  and  (Cb  -  9)  NLaR  %0-x 


The  integrals  in  Equation  (Cb  -  8)  and  (Cb  -  9)  are  of  the  f  crm: 


J  ->-< *>*)  "§■<>» 

(IV  -  1) 

The  integral  may  have  been  evaluated  to  such  a  value  near  x  that  for  the 
remaining  interval  from  %•  to  x  it  can  be  replaced  by: 

J  ■  <§■>.  A  -  (£>‘]b «» 

(IV  -  2) 

with  the  average  temperature  gradient  (dT^/d  $  )  in  this  interval  determinable 
with  sufficient  accuracy  by  an  estimate. 

fy  a  change  of  the  variable  %  to:  y  *  x-  §  and  introduction  of:  e  ■  x-)a 
obtain: 


J  -<§■). /[1  -  u -(?*] b  v 

9  O 

(IV  -  3) 

Expanding  (l-(y/x))a  into  a  series  and  neglecting  terms  higher  than  first  order 
gives: 


(1  1 


Introduction  into(l?  -  3)  and  integration  finally  results  in: 


(IV  -  4) 


HT  ab  a6*1 

1  -  <$■>.  s5i  V 


(IV  -  5) 


we 
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The  location  of  $«  can  always  be  chosen  such  that  the  approximation  (ZV  -  4)  is 
fulfilled  with  sufficient  accuracy. 

ThenCtV  -  5)  gives  the  contribution  of  the  interval  between  and  x  to  the 
integrals  (Cb  -  8)  and  (Cb  -  9)* 
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Figure  2.  Calculated  Laainar  Flat  Plate  Recovery  Factor  as  Function  of  Pr  Number 
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Figure  3.  Properties  of  Air  as  Function  of  Temperature 
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Figure  4.  Calculated  Skin  Friction  Coefficient  for  Laminar  Boundary  Layers  on  Flat  Plate. 


Figure  5 


Temperature  Profiles  and  Reference  Temperature  in  High 
Speed  Boundary  Layers. 
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Figure  6.  Calculated  Laminar  Recovery  Factors  as  Function  of  Mach  'lumber 
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Calculated  Laminar  Fnthalpy  Recovery  Factors  for  Flat  Plate 
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Figure  3.  Pwatio  of  Friction  Factor  to  Fnthalphy  Stanton 
"umber  for  Lminar  Flow  over  Flat  Plate 
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Figure  9.  Ratio  of  Actual  to  Constant  Property  Turbulent  Friction  Factor  for  Flat  Plates 
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Figure  10.  Calculated  Turbulent  Friction  Factors  for  Flat  Plates 


Figure  12.  Laminar  H°at  Transfer  on  Flat  Plate 
with  Variable  Wall  Temperature. 


84 


X  <— •* 

5  — 


Figure  13.  Calculation  of  Heat  Transfer  Along  Flat 
Plate  with  Variable  Wall  Temperature 
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Figure  14.  Flow  Regimes  in  Kach-Reynolds  Field. 
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Figure  15.  Stanton  Numbers  for  Heat  Transfer  in  Free  Molecular  Flow  Regime 
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Figure  16.  Temperature  R°cov°ry  Factor  for  Free  Folecular  Flow  Regime. 


Figure  17.  Average  Skin  Friction  Factor  for  Flat  Plate  in  Rarified  Air 
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Figure  18.  Average  Stanton  dumber  for  Heat  Transfer  on  Cones  in  Rarified  Air. 
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Figure  19.  Maximum  Temperature  T^aoi  Recovery  Temperature  Tr  in  Bouniary  Layer  on  Flat  Plates. 
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20.  Reynolds  ,Tunber  fo 
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Figure  22.  Calculation  of  Reference  Temperature  for  Laminar  Flow  Along  Flat  Plate 
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